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Abstract 



The review is devoted to the theory of collective and local pinning effects in various disordered 
non-linear driven systems. Common feature of both approaches is the emergence of metastability. 
Although the emphasis is put on charge and spin density waves and magnetic domain walls, the 
theory has also applications to flux lines and lattices thereof, dislocation lines, adsorbed mono-layers 
and related systems. In the first part of the article we focus on the theory of the it collective pinning 
which includes the equilibrium properties of elastic systems with frozen-in disorder as well as the 
features close to the dynamic depinning transition enforced by an external driving force. At zero 
temperature and for adiabatic changes of the force, the dynamic depinning transition is continuous, 
the correlation length is large, the behavior is dominated by scaling laws with non-trivial static and 
dynamical critical indices. The application of functional renormalization group methods allows for 
the detailed description of both equilibrium as well as non-equilibrium properties. The depinning 
transition is also characterized by the appearance of it new scaling laws. Thermal fiuctuations smear 
out this transition and allow for a creep motion of the elastic objects even at small forces. The 
application of an ac-driving force also destroys the sharp transition which is replaced by a velocity 
hysteresis. 

The second part of the review is devoted to the picture of the local pinning and its applications. 
Local theories apply in the region where correlation effects are less important, i.e. not too close 
to the depinning transition, at low temperatures, at high enough frequencies or velocities. The 
inclusion of plastic deformations results in a rich cross-over behavior of the force-velocity relation 
as well as of the frequency dependence of the dynamic response. Being easily affected at higher 
frequencies or velocities, the local pinning becomes the easily accessed source of dispersion, relaxation 
and dissipation. The picture of the local pinning can be effectively used to explain experimental data: 
qualitatively and even quantitatively. The advantages come from the explicit treatment of metastable 
states, their creation and relaxation, their relation to plasticity and topological defects. The local 
pinning recovers and exploits new elements of the energy landscape such as termination points of 
some branches or irreversibility of other ones related to generation of topological defects in the course 
of sliding. It provides also a clue to quantum effects describing the quantum creep as a tunnelling 
between retarded and advanced configurations. 
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1 Introduction 



Many ordering phenomena in solids are connected with the emergence of a modulated structure. Ex- 
amples are charge or spin density waves PEIEIEEIEIIZIIH], Wigner crystals mEl, flux hue lattices 
|l()l , incommensurate phases of adsorbed monolayers ^2] etc. These structures often interact with 
imperfections frozen in the solid leading to pinning phenomena which change drastically the statics and 
dynamics of the modulated structure. 

In other cases topological defects like isolated flux lines in supcrfluids or superconductors, dislocation 
lines in solids, domain walls in magnets etc. appear as a result of competing interactions, external 
fields, or in the process of fabricating the material. Pinning of the motion of these objects is often 
required if one wants to exploit physical properties of the condensed structure. The pinning of flux lines 
prevents dissipation from their frictional motion in superconductors, pinning of dislocations prevents 
plastic deformations of a solid etc. [131 ITTl I15| . The generality of this approach has been anticipated 
already long time ago but strong similarities between different systems mentioned above have been 
worked out in detail only later on. 

The goal of the present article is to give a presentation of unifying concepts in the theory of pinning 
phenomena. To make the ideas more clear, we will not go too much into the details of specific systems 
but we shall stress the generality of the approach. Some aspects of pinning have been considered in the 
past in great detail for type-II superconductors and we refer the reader to these articles for more details 

nniin]. 

A perfectly rigid object of any dimension, e.g. a straight flux line, a planar domain wall or an 
undistorted charge density wave will never be pinned - the pinning forces acting on different parts of the 
object cancel each other. More accurately: the resulting total pinning force is of the order of the square 
root of the volume of the object. It is therefore necessary to consider the distortions of these objects 
under the influence of pinning centers which are often of elastic type. We will therefore speak about 
"elastic systems" if we refer to arguments which are essentially correct for all systems mentioned above. 

The delicate question of pinning by randomness stays already for almost three decades since the earliest 
proposals of Larkin ^21, through theories of vortex lattices JHl and density waves ^112011^. The last 
decade has brought new insights, largely provoked by the studies of vortices in high-Tc superconductors 
|1U[ . In particular new understanding was reached detecting the possibility of a quasi long-range 
ordered (Bragg) glass phase with an algebraic decay of structural correlations in a disordered system 
|22l 1231 El] and by the understanding of the role of mctastablc states in non-stationary effects. Essential 
progress was due to new advanced methods like the functional renormalization group method |25l I2f>j for 
the collective pinning problem. Still full details are not completely accessible and applications are not 
always straightforward. The picture of the local or weak pinning was also revised and developed through 
the last decade in conjunction to plasticity and a role of topological defects. 

The undisputable domain of the collective pinning approach are the effects occurring on large length 
and time scales. These are recovered e.g. in studies of long time evolution, the low frequency response 
to external forces, the creep below the zero temperature depinning threshold, and the region around the 
threshold field. The collective or weak pinning forces come from elastic interference of many impurities 
|18[l21ll^ll9| . Their characteristic features are: large correlation volumes - of the order or beyond the 
Larkin-length, high energy barriers between mctastablc states, huge relaxation times and small pinning 
forces. The collective pinning determines the threshold fc of the driving force to initiate the sliding quite 
similar to conventional rest friction |27j . The comprehensive picture of the collective pinning will be given 
in the first part of this review in Sections 2-5. 

Complementary insight can be obtained within a frame of the so called strong or local pinning which 
manipulates with only a few of metastable states. This will be the topic of the second part of this review 
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(Sections 6-12). This simple but transparent approach offers some effects which have not been noticed 
yet - and still are not fully accessible - within the very complicated picture of the collective pinning. The 
local or strong pinning comes from rare metastable pinning centers which provide finite barriers, hence 
reachable relaxation times. Being easily affected at higher frequencies or velocities, the local pinning 
becomes the easily accessed source of dispersion, relaxation and dissipation. This part will summarize 
the local pinning approach to time dependent properties of sliding superstructures. We shall follow the 
scheme of |2H1 EHl EOl ED for the theory of the pinning-induced metastable plastic deformations partly 
due to creation of dislocation loops or lines. Depending on the pinning potential strength, we shall find 
several regimes of the behavior of local deformations in the course of a displacement with respect to 
the impurity site. The key observation is that the local state at the pinning center - the impurity - 
may be either unique or bistable. The bi-stability can be either restricted or unrestricted (i.e. preserved 
throughout the whole period of sliding). In the latter case this leads to the generation of diverging pairs 
of dislocation loops - or 2tt~ solitons in a quasi 1-dimensional picture. On this basis we can obtain 
contributions to the pinning force resulting in the sliding velocity-driving force characteristics and in the 
frequency dependent response. 

Special applications will be devoted to the best studied examples of domain walls and CDWs where we 
shall particularly address two commonly observed experimental features: the totally nonlinear current- 
voltage curve and the anomalous low frequency low temperature behavior of the dielectric susceptibility. 

2 Equilibrium properties of elastic objects in random environ- 
ments 

In this section we summarize the equilibrium properties of elastic objects in a random environment. 
There are two main classes of systems: The first one refers to non-periodic objects like isolated domain 
walls, flux or dislocation lines. These can be generalized to the so-called elastic manifolds models. The 
second class of models is periodic, like charge density waves, flux line lattices or Wigner crystals. These 
are subsumed under the expression periodic media. 

To make the notation simple, we will describe the distortions of the elastic systems from perfect 
order by a scalar displacement field w(x), its generalization to A^- vector or more complicated fields is 
straightforward. Correspondingly we will mostly use the terminology of domain walls and charge density 
waves, respectively. Where peculiarities in systems with N > 1 may occur we will mention them however. 



We consider a Z)-dimensional elastic object embedded into a host medium of a dimension d. Since the 
medium includes quenched disorder, the energy can be written in the form 



Here D denotes the internal dimension of the object {D = 1 for isolated flux or dislocation lines, D = d—1 
for isolated domain walls, D = d for charge density waves, flux line lattices and Wigner crystals etc.). 
Note, that D is in general different from the space dimension d {D < d). (The dimension TV denotes the 
number of components of the displacement field u. For isolated flux lines or domain walls N = d — D 
whereas for fiux line lattices N = 2 and iV = 1 for CDWs.) C is an elastic constant and in general a 
tensor, e.g., in flux line lattices at least three elastic constants are necessary for the description. In some 
cases the elasticity is non-local on intermediate length scales - as for flux line lattices on scales smaller 



2.1 Models 




(1) 
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then the London penetration length - or even at all length scales as for dislocation lines. The elastic 
constants will also show a temperature dependence since they have to vanish at the transition where the 
structure melts, e.g. at the Peierls transition for CDWs. / denotes an external force density which acts 
on the object. 

The random potential Vr(x, u) results from the coupling of the elastic object to the impurity potential 
u/j(x, z), which is generated by the frozen impurities, fluctuating exchange constants etc.: 

Vr(x,u)= J (i'^"^zui?(x,z)p(x,z,u) , (2) 

/9(x, z,u) represents the density of the elastic object and will be specified below, x = (.ti, ...xd) and 
z = {xd+1, ■■■,Xd) denote a D- and {d — Z?)-dimensional position vectors parallel and perpendicular to 
the main orientation of the object (e.g. the domain wall or the flux line), respectively. For density waves 
and flux line lattices D = d and hence there is no perpendicular coordinate z. 

The average over disorder - which replaces the average over the (infinite) sample - will be denoted 
{...)r. Without restricting the generality we will assume that (u_r(x, z))ij = - with other words, we 
incorporate the effect of the averaged disorder potential into the bare parameters of our model and 
consider only its fluctuations. A less general but still reasonable choice is that Uii;(x, z) is Gaussian 
distributed and short range correlated (with a correlation length /) such that it is characterized by its 
second moment: 

(wH(xi,Zi)wij,(x2,Z2))_R = w|j(5(xi - X2)5(zi - Z2). (3) 

If we rewrite the random potential ^^^^(x, z) as a sum of impurity potentials 

vr{x,z) = ^ ViS{x - Xj)(5(z - z.,) - ViUimp (4) 

i=l 

then the disorder average is defined as 

iV^P . ,D Jd-D r 

{■■■)n- n J V J ^^'^(^^)--- (5) 

-00 

Here Nimp, V and riimp denote the total number of impurities, the volume of the sample and the impurity 
concentration, respectively. p{Vi) is the normalized probability distribution of the potential strength Vi 
with = / dl/, . . p{V{) and = V^n.mp. 

The correlations of the random potential vr are always short ranged provided the impurities are short 
range correlated. (5(x), (5(z) are J-functions which may be smeared out over a length scale of the order I. 
Since we will use it later we introduce here also the correlation function 

(Vr(xi,Ui) VR(x2,M2))fl = i?(ui,M2)(S(xi - X2). (6) 

Because of the gaussian nature of the disorder correlations the disorder averaged physical quantities like 
the free energy depend only on the correlation function i?(ui, U2). The relevant contributions of i?(ui, U2) 
are those which depend only on the difference ui — U2 on which we will concentrate in the following. For 
|u| <^ I, R[u) is quadratic in u (with i?„„(0) < 0), for larger u different cases have to be distinguished 
(see below). 

In the following we consider some specific examples: 



6 



(i) Domain walls in magnets : Here u(x) describes the displacement of the domain wall from a planar 
reference configuration, D ~ d—1 and = 1. The stiffness constant C is finite only above the roughening 
transition temperature Tr, for T < Tr the elastic description breaks down in the absence of disorder. 
But disorder leads always to roughening of the interface, even at T = |321 El |^ . The driving force 
density is directly related to the magnetic field B hy f = 2^bB. If the random potential results from 
fluctuations of exchange coupling between the spins (this is the so-called random bond case) , which couple 
only to the domain wall, then /9(x, z, u) ^ d{z — u(x)). In some cases the width of the domain wall may be 
large compared to the lattice spacing and the (5-function has to be replaced correspondingly by a smeared 
out profile function. The correlations of Vr(x, u) are then short range in x and have a correlation length 
corresponding to the maximum of / and the domain wall width in the it-direction. To keep the notation 
simple we will denote this maximum in the following also by I. R{u) is then a 5-function of finite width 
of the order ~ I. The extension of this model to arbitrary D and N (i.e. D + N is not longer equal to d) 
is called the random manifold model. 

If disorder comes from random field impurities then the domain wall couples to the disorder in the 
domains. In this case p(x, z, u) ^ 0(u(x) — z) which gives the Zeeman energy with a constant magneti- 
zation in each domain, w^Xx, z) represents now an uncorrelated random magnetic field. In this case of a 
non-local coupling to the disorder it can be shown that R{u) « — |it| for |u| ^ I |35j . 



(ii) Isolated magnetic flux lines or dislocation lines: u(x) — > u(x) denotes now a N = d—1 component 
displacement field and D ^ 1. The random potential yR(x, u) is short range correlated both in x and 
u and hence i?(u) is again a smeared out 5-function. For flux lines, / is given by the Lorentz force 
f = ij X b<I>Oj where j denotes the transport current. $0 = hc/2e is the flux quantum and b the local 
direction of the magnetic field. In the case of dislocation and vortex lines the elastic energy is non- 
local In particular, for dislocation lines one finds after Fourier transformation a weakly momentum 
dependent elastic modulus C{k) = —C\og{aok). oq is of the order of the lattice constant. The force 
acting on the dislocation line is the Peach-Kohler force |3()| . 

(iii) Charge density waves: The condensed charge density can be written as 

p(x,(/j) =po(l + Q"'V</j) + Pi cos (Qx + (^(x)) + ... (7) 

where Q denotes the wave vector of the charge density wave modulation and the dots stand for higher 
harmonics. The first factor describes the density change due to an applied strain. The phase field 'p(x) 
is related to a displacement u(x) of the maximum of the density by 

(p(x) = -u(x)Q, (8) 

which we will use from now on as the relation between u and ip. 

Since charge density waves carry an electric charge, the stiffness constant C shows in general a strong 
dispersion due to the long range Coulomb interaction 

Ck' ^ C|| kj + C^kl + Cd^^h/^ , (9) 

which we will ignore except for special applications. Here A denotes the screening length. If A diverges, 
the system may be considered effectively as a 4-dimensional one with k\\/k playing the role of a fourth 
dimension |?7| . 

Corresponding to the two (^-dependent contributions proportional to po a-nd pi in |(7J) there are two 
contributions in ^ which sometimes are referred to as forward and backward scattering, respectively. 
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The resulting correlator R{u) is periodic in u = —ip/Q with the periodicity 2n/Q. Indeed, for charge 
density waves we obtain from equations @ and 

Ri^)'x^^pl{Q-'V^f + ^pI cos if. (10) 

Here we have neglected strongly oscillating terms (which average to zero) as well as terms which can be 
included into the elastic energy. Finally, the external force / is given by the applied electric field, f ^ E. 

For a more detailed discussion of other systems like flux-line lattices, Wigner crystals etc. we refer 
the reader to the appropriate literature |1UII11[I551I39[ I4U) . 



2.2 Basic properties of disordered systems 

To get preliminary information about the influence of disorder on the elastic object and the relevant 
length scales we consider first small distortions w(x) around the state of perfect order ?i(x) = 0. The 
energy can then be written as a series expansion in it(x): 

j d'^xi^civuf + Vr{^,q) + Vr.u{^,q)u + ---Y Vr^u^^. (11) 

Here we assumed that the distortions are small such that we can neglect for the moment higher order 
terms in the expansion of Vr{'x., u). 

The necessary condition for the ground state of the Hamiltonian follows from the variation of Ti with 
respect to u 

^ = -CV^t + FH^„(x,0) = 0, (12) 

which is the Poisson equation known from electrostatics [u and VR_ti(x, 0)/C playing the role of the 
electrostatic potential and the charge distribution, respectively). 
Its solution is given by 



^j,, Fi^,„(x',0) 



From ((T^ we obtain with the help of ® 



where i?M«(0) denotes the second derivative of the correlator at u = 0. At this point we want to mention 
already one problem with the use of the approximation H12() : since the random forces Vr.„(x, 0) do not 
depend on u, any constant can be added to the solution H13|l . In other words, the object could be moved 
through the random environment without any change in energy. Thus the linear approximation misses 
barriers and metastable states. 

As can be seen from (I14f) . the relative displacement in two different points 0,x increases with their 
separation |x| below the critical dimension Dc = 4, the elastic object is called to be rough. We can rewrite 
H14|l in a more general form as 

wr{L)^i(^^^\ (15) 

where we introduce the roughness exponent C^. In our present calculation Q = ^rf = (4 — -D)/2, where 
the subscript RF stands for random force corresponding the expansion of Vr(x, u) up to a force term in 
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equation Ijllf) . The characteristic length scale Lp 

is called the Larkin-length |17j in the context of flux line lattices, the Fukuyama-Lee length |20| in the 
context of charge density waves, the Imry-Ma length for random magnets etc.. On this length scale 
the displacement is of the order of the correlation length I of the random potential, i.e. the displacement 
field can choose between different energy minima and hence meta- stability appears. For weak disorder, Lp 

is large compared with the mean impurity distance n^^^ and hence pinning phenomena arise from the 
collective action of many impurities. In the situation of strong pinning, which we will discuss further in the 
second half of this article, individual pinning centers lead to strong distortions already on the scale shorter 
than the distance between impurities n^^^'^. In the next chapters we will always assume Lp ^ n^^p . 
The Larkin-length will then play the role of an effective small-scale cut-off for the phenomena considered 
in the following. 

Unfortunately, the result 1)14(1 is not applicable on length scales much larger than Lp since equation 
p2|l is only a necessary condition for the ground state: it is the condition for a saddle point, not only for 
the absolute minimum. It is generally believed that the elastic object has a unique (rough) ground state 
|42|. But for distortions u > I the system has in general many local minima and perturbative methods 
break down. A nice demonstration of the breakdown of perturbation theory in systems with several 
energy minima has been given by Villain and Semeria |43| . In this case more elaborate methods like the 
renormalization group approach have to be applied. This exceeds the scope of the present short review, 
but a schematic presentation of the method is given below. 

The result of this approach is that the roughness w{L) can still be written in the form (|15() but with 
(^RF replaced by a non-trivial roughness exponent C, with < ^ < 1. (For C > 1 | Vu |~ L^^^ diverges on 
large length scales and the elastic approach adopted here breaks down.) 

For periodic media, in general, effects of disorder are weaker. In the case of elastic manifolds it 
may pay off for the elastic object to make a large excursion to find a larger fluctuation in the impurity 
concentration, which lowers the free energy. On the contrary, for a periodic medium, the distortion has 
to be at most of the order of the period of the medium to reach a favorable interaction between the 
impurity and the medium. Hence we expect that periodic systems belong to another universality class 
with different, smaller exponents. The actual result happened to be more drastic: the roughness exponent 
( for periodic media is zero corresponding to a logarithmic increase of wr{L) |22I \2'6\ 1^ 144) (for further 
details see Section IT^ . 



The average (...)t over thermal fluctuations does not change the asymptotic behavior of the correlation 
function 

«;|j(|x|) = ((«(x)-«(0))|)^(x|xpC (17) 

which vanishes apparently in the absence of disorder. It is this correlation function which dominates 
the structural properties of the system. But the temperature will have a drastic effect upon the time 
dependent properties. 

Since we want to consider later on the motion of the elastic object under an external force / it is 
instructive to characterize the energy landscape the elastic object is exposed to. As a first step we consider 
changes in the free energy Tr (which depends on a particular configuration of the disorder) by going 
over to a different disorder configuration. Such new disorder configuration may be created in some cases 
by applying external forces (compare with equation |(S7|)) or changing the boundary conditions. Sample 
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to sample fluctuations of the free energy of a region of linear extension L arc then expected to show also 
scaling behavior described by a new exponent x 

i^[TR{L)~{TR{L))Rf)'^^ ^T,(^^y = F{L), Tp^CPL^-^ x^D-2 + 2Q. (18) 

where Tp denotes a characteristic energy scale see e.g. ^Hl^lEIEHl- The scaling relation between x 
and C can be understood from the fact that in the free energy the elastic and the random part of the 
energy have to be of the same order and therefore it is plausible that the scaling behavior of !Fr{L) 
can be read off from the scaling of the clastic energy, which scales with the exponent x- An illustrative 
one-dimensional example of the free energy fluctuations is considered in Appendix 1X1 

Further information about the system comes from a second correlation function uiydxl) which de- 
scribes the response to a local force which couples to (m(x) — m(0)): 

«;|(|x|) EEE (([u(x) - u{Q)f)T - {u{^) - u{Q))l)^ ^ T|x|2-^/C + const. (19) 

Any pair correlation function of it(x) can be expressed by a combination of Wfl(|x|) and wt(|x|). Equation 
l)19|l is (exactly) the same result as in a non-random system and a consequence of a statistical "tilt- 
symmetry" of the system ^^l- Apparently, u;t(|x|) vanishes at T = 0. The result p9|l can be related to 
the static susceptibility 

X=(^WO))t^^ = ^ / d''x{{u{^)umT^{u{^))T{umT)R (20) 

where / is the force conjugate to u. 

To discuss the result for wt(|x|) further we consider the case when u(0) is fixed to zero such that 
uiydxl) = (((m(x) - (M(x))T)^)y)^ = (((5u^(x))^)^. The quantity (5it(x) = m(x)- (u(x))t describes the 
fluctuations of m(x) around its thermal average in a given random environment. Naively one could assume 
that these fluctuations are restricted to a narrow valley along the ground state such that the fluctuations 
should not increase with |x| = L. The result H19() however suggests a different picture: besides of the 
ground state there are rare excited states which are very different in configuration from the ground state, 
(5w(x) ^ 1{L/Lp)'^, the energy of which differs only by an amomit of the order AE < T from the ground 
state (compare with Figure^). 




X 



u 



Figure 1: Ground and low energy excited state of a ID domain wall 

These excited states could become indeed true ground states if we would change the random potential 
locally (i.e. in the neighborhood of the initial ground state) in an appropriate manner. Since at T = 
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the free energy is given by the ground state energy we have to expect that rare excited states differ in 
energy as the sample to sample variations of the free energy. 

To calculate the fluctuations of Su we make therefore for the probability distribution of the energy of 
the excited states of an elastic object of linear extension L the scaling Ansatz 

PiAE, L) = F-\L)p{AE/F{L)) (21) 

where F(L) denotes the sample to sample variations of the free energy (compare with equation H18|l 'l and 
p{x) is an unknown normalized function with p(0) = 0(1) > 0. With this Ansatz we find for the average 
fluctuations at low temperatures T ^ F{L) |47l I48j 

{\Su{^)nT,R « ^" [y-J I d{AE) P{AE, L)-r- i^-J (22) 

which gives for n = 2 the exponent — % = 2 — D appearing in (|19|) . 

The success of this approach gives us the possibility to calculate also the specific heat 

r) 8 f°° /"^ 

c{T) = L^^ — {AE)a^— dLv{L) d{AE)-AE-P{AE,L). (23) 

•J Lp •J 

Here Lq and i^iL) denote the system size and the size distribution of the rare low energy excited states 
on the scale L, respectively. In writing down (|23|l we have decomposed the system in (Lq/L)^ blocks of 
linear extension L. Each of them gives a contribution to the internal energy of the order of the integral 
on the r.h.s. of 1)23(1 . Then we have to sum over the contributions from all length scales L, Lp < L < Lq. 
The smallest scale is clearly given by the Larkin scale Lp. The next independent contributions come 
from excitations on larger scales bLp,b^Lp, ...,b"Lp etc., 6 > 1. (5 has to be chosen in such a way that 
excited states on scale b'^'^^Lp cannot be reached already on scale b'^Lp, which requires b'' > 2.) The 
total number of scales is given by ln(Lo/ip)/ln6. If we replace this sum over n by an integral over L, 
we obtain a factor (the integration measure ) di(Lln&)-i(^)^ and hence iy{L) cx j^]]. (Roughly 

speaking we could say that we integrate over all momenta d'^k with k ^ L^^.) 

In principle, the distribution P{AE, L) may also depend on the temperature. If we ignore this 
unknown dependence, we get 

since the integral is dominated by small L and p{T/Tp) k, p(0). This approach gives at T <C Tp a specific 
heat linear in T which is similar to physics of common two- level systems in amorphous solids |51| . This 
analogy also builds a bridge to the local pinning picture which thermodynamically is equivalent to the 
case of |51j . Within the collective picture, similar results have been obtained, in another way, recently in 
|52j . The precise value of b remains unknown in this approach but has no influence on the temperature 
dependence. Measurements of the specific heat on a finite time scale t will lead to a reduced value of 
c(T) since not all local energy minima can be reached by thermally activated hopping. We will give a 
time dependent correction factor (which takes this fact into account) at the end of Section 3.2. 

So far we have considered mainly the case of weak pinning. The large scale properties (L > Lp) 
are dominated in this case by density fluctuations of the disorder where many impurities are involved. 
We want to stress that in the case of strong pinning on very large scales pinning phenomena are again 
dominated by density fluctuations of the impurities. Strong pinning on small scales is however very 
different from weak pinning and will be further discussed in the second part of this article. A one- 
dimensional CDW-model with strong pinning ^1|53] is considered in Appendix IbI 
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2.3 Results from the renormalization group method 

In the rest of this section we give a short account of the apphcation of the renormalization group (RNG) 
[5^ ESI approach to elastic systems in random environments EZ] ■ The RNG method starts from 
a Fourier decomposition of the displacement field u(x) = X^k '"kS"*''''- The first step of the procedure 
consists in the elimination of the short wave length degrees of freedom u^. = with k in the momentum 
shell Ao/6 < |k| < Aq, 6 > 1, from the partition function. Aq = 27r/ao denotes a microscopic momentum 
cut-off. The result of taking the trace over can be written again in the form of a Boltzmann factor with 
a new effective Hamiltonian for the remaining degrees of freedom. If thermal fluctuations are irrelevant, 
as for D > 2-dimensional systems, this procedure reduces to the problem of finding the values ~ ii^ 
which minimize the energy, keeping all Uj, = with |k| < Ao/6 fixed. Plugging these values into 
the Hamiltonian we get the new effective Hamiltonian which contains less degrees of freedom Since 
the scale baQ on which the components describe displacements is small, we can expect that there is 
only one minimum as a solution for such that the application of perturbative methods is allowed. 
To avoid further misunderstanding we stress here that we always assume in this Section that the system 
may reach thermal equilibrium, even if the time scales are huge. Besides of corrections to terms already 
present in the initial Hamiltonian also new terms may be generated in this procedure, their precise form 
depends on yR(x, u). The concrete implementation of this procedure is in general difficult and rests often 
on approximations valid close to certain critical dimensions. 

The second step of the renormalization group procedure consists in a rescaling of length, time and 
fields according to 

x = x'6, t^t'b^, u(x) = u'(x')6'^ (25) 

with so far unspecified dynamical exponent z (> 0) and roughness exponent {> 0). After the first step 
of the RG-transformation the remaining minimal length scale was Axmin ~ aob which, after rescaling 
according to goes over into the original minimal length Ax'^^^ = Axniin/6 = ag. 

Using rescaling (|25|) in equation (Q, the elastic energy in the new coordinates obtains a factor 6^^^+^^ . 
Since in statistical physics the Hamiltonian appears always in the combination 7i/T, we can absorb this 
factor in a rescaled temperature 

T' = r62--D-2C (26) 

The rescaling has also to be applied to the second term in Having introduced T', Vr(x, u) is replaced 
by 62-2Cv^(6x',m'(x')5«) = V^(x',u'). The correlator of is then given by 

(FA(x;, z.;)F^(x^, 4))^ = &4^4f-D5(^, _ )^ (^C(„; _ ^^)) . (27) 

As mentioned already, equations (|26|l and l|27|) do not represent the whole change of T, R under the 
renormalization, since the new terms produced in the first - non-trivial - step of the procedure will 
generate contributions to R' and T' . In general the whole function R(u) is transformed in a non-trivial 
way. We will not discuss here the derivation of this transformation but only present the result. The 
interested reader is referred to the original articles 1561 157] . For infinitesimal changes oi b = 1 + Sb 

one obtains a continuous flow of these quantities. To lowest order in 4 — D one finds with 6b ~ 6L/L: 

(4 - D - 40R{u) + CuRuiu) + l-Ruuiuf - Ruu{u)Ruu(.0), 

^ (28) 

{2~ D~ 2C)T + ... 

where i?„(w) = -^R{u) etc. To make the notation easier we have adsorbed here the coefficients of the 
quadratic terms in R into the redefinition of R{u). 



dRju) 

dh-iL 
dT 

d\nL 
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A word of caution seems to be indicated at this point: the rescahng (|25() is a matter of eonvenience, 
not a physical necessity. If one does not rescale at ah, the first step of the renormahzation group still 
keeps the whole information about the large scale behavior of the system. In this case one obtains the 
effective physical quantities (those which one observes in experiments) on the corresponding length scale 
L. 

For systems with uncorrelated disorder (Vr(x,u) is a random function of both arguments) there is 
an important simplification since there is no renormahzation of T . Indeed, the full stiffness constant 
C can be measured by replacing the periodic boundary conditions used so far by u[L,X2, ■■■,xd) = 
u{Q,X2, ...^xd) + Su, that is applying an overall strain ^ in xi^direction. The effective elastic constant 
Ceff follows then from the free energy J-b.{L, u) 

Ceff^ L-"+^d^TR{L,u)/d5u^ . (29) 

The change of the boundary conditions can be compensated by introducing a new variable u{yi) — 
w(x) — 5uxi/L. Changing from m to u in the Hamiltonian adds to the elastic energy a constant 
contribution SH. = ^L^^^Su^ and changes the random potential into / d'''xVfj{x,u{x) +Suxi/L). Since 
Vr(x, u) is by definition a random function of its arguments, V/{(x, m(x) + 6uxi/L) can be replaced 
by a random potential ^(x, u(x)) with the same statistical properties. Thus, the only change in the 
Hamiltonian is given by the constant STi, = 5J-r{L, u) from which one concludes with H29|) Ceff — C . 
Since T appears only in the combination C /T , there is also no renormahzation of T . 

D i 



one flat phase: ^=0 



4 




Figure 2: The D — N plane characterizing different elastic manifolds in a random medium (a ~ N). For 
D > A weak disorder and for D > 2 thermal fiuctuations are irrelevant. For A > D > Dn (i.e. Xb. > 0) 
weak disorder is relevant leading to a non-zero roughness exponent. For Dn > D a, thermal depinning 
transition exists: at low temperature weak disorder leads to a non-trivial roughness exponent C whereas 
in the high temperature phase Ct = (2 — D)/2. 

The flow of R{u) and T as given by equation (|28|l will terminate in stable fixed points (other forms 
of flow like limit cycles are excluded by general reasons), which characterize physical phases. To flnd a 
finite fixed point of the disorder (if it exists) we have to tune the value of the exponent C in order to make 
elastic and disorder energy scale in the same way. The flxed point condition thus delivers the value of 
the exponent C,. Different flxed points have in general different values of C,. If several fixed points exist 
- this is the typical situation - the domain of attraction of a fixed point determines the phase boundary. 
Clearly, the flow close to a fixed point depends on the value of C, at this fixed point. 
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D>4 



2<D<4 




T 



Figure 3: Schematic phase diagrams of impure elastic systems as a function of space dimension for 
(left) non-periodic media (single flux lines, domain walls) and (right) periodic media (flux line lattices, 
charge-density waves, Wigner crystals), vji and T denote the strength of the disorder and temperature, 
respectively. Fixed points characterize the properties of a phase, their domain of attraction ends at the 
phase boundary. For realistic phase diagrams further degrees of freedom may become relevant (disloca- 
tions etc.). In Z? > 4 dimensions besides of the fixed point v*^ = there may be another strong disorder 
fixed point. 



If there is no disorder, i? = 0, the roughness exponent is C = Ct = (for -D < 2) according to (|28|l 
and each value of T is a fixed point. Since we exclude negative roughness exponents (there is always a 
constant contribution to the roughness), we put Ct = for D > 2 and hence T is transformed to zero. 
Switching on the disorder only the linear terms in R in equation H28|l matter as long as the disorder 
remains weak. If the initial function Riu) shows a simple power law behavior, uRu{u) oc —aRiu) 
[a = —1 for random field and a = N for random manifold systems, respectively), we get d\nR/d\nL — 
4— D — C(4-ha) = Afl,(C). Using as the initial value for Q its value at the thermal fixed point, Q = ^t, wc see 
that the disorder grows provided 4— {A + a)(T — D > 0. To get a new stable fixed point wc have to choose 
then a new value for ^. In this situation we have to take into account also the nonlinear contributions in 
equation (|28|l . The resulting fixed point function R*(u), which determines also the value of C, depends 
on the initial function of R{u) and can be found often only numerically |25[ I26| . A characteristic feature 
of the R* ( u) is the cusp singularity, i?*^ (u) — (0) oc | u | for small | u | , which appears on length scales 
L > Lp and is related to the appearance of metastability (see below) . 

The true roughness exponent for manifolds in random bond systems can then be written as an expan- 
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sion in e = 4-1). For interfaces in random bond systems one finds ( = 0.2083e + 0.0069e2 + O(e^) EElEHl- 
It has been suggested that Cfls can be written in closed form as = 4+n^{d n) where 1) = 0.5, 
i/(l,2) = I and 1^(1,3) = ^ 

For random field system (rf = for 2 < £> < 4 [201 • 

For Xr{Ct) > 0, T/C is always transformed to zero, thermal fluctuations are irrelevant. In the 
opposite case Xii{C,T) < weak disorder is irrelevant. Increasing the disorder strength, there is a phase 
transition between a disorder dominated phase at low temperatures and a high temperature phase where 
thermal fluctuations wipe out the disorder. The relation Xr{Ct) = define the dimension D{N) in 
FigureEl 

The discussion applied so far to the so-called random manifold models, which are given by a non- 
periodic disorder correlators R{u). For periodic media the correlator R{u) is periodic with period b and 

this applies also to the fixed point function R*{u) ^i^R^if) (I^t) )• However the term (uRu{u) 



36 ^ ^ ^ 

in equation H28() violates the periodicity, from which one has to conclude that for periodic media = 
corresponding to a logarithmic increase of the roughness |22l 1231 1241 1441 IM] : 

wl{L)^i4-D)lHn(^^y (30) 

For some applications, e.g. for the calculation of the correlation function wr{1), (but not for studying 
the effect of energy barriers) one can derive the results for L > Lp from the random force model ()12|l but 
with a modified correlator for the random forces Vr^„(x) 

{VnA^,0)VnA^\0))n = J rf^fce"-(--')^--^ -|°],^^^_^ (31) 

where i?*^ (x 4 — D. As can be seen from (|31|l . the force correlations on length scales L ^ k^^ <C Lp 
behaves as i?„„(0)(5(x — x') as in the random force approximation of Larkin [17], whereas on length scales 
L ~ ^ Lp there is a long range contribution decaying as |x — x'|^^. 

Using the result equation (|30|) in the structure factor S'(k) one obtains a smeared (diffuse) Bragg-peaks 
of finite width, 5(k) ~ |k — G|~'^+'"=, despite of the fact that the system is dominated by the influence 
of disorder, hence the name Bragg-glass has been coined |22l 1231 I61| . 

In D > A weak disorder flows to zero, i.e. it is irrelevant, but for sufficiently strong disorder a separate 
strong disorder fixed point for R may exist. It should be noted that periodic media allow the existence 
of topological defects like vortices or dislocations which may destroy the Bragg-glass phase. However it 
can be shown that for weak enough disorder this phase survives in d = 3 dimensions [5^ ICT IMI ICT ISS) . 

In D = 2 the T-axis (i? = 0) is a line of fixed points corresponding to C = 0. For random manifolds 
disorder is always relevant, i.e. C > ^nd there is a non-trivial fixed point at T* = 0, R* > (see the 
left side of Figure |2J). For periodic media (like flux line lattices in a thin film or Wigner crystals) the 
situation is more complicated: there are two phases both with C = separated by a phase transition at 
T = Tg (see the right side of Figure EJ. In the high-temperature phase 

w'iL)^^\n(^^y (32) 
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where as in the low-tempcraturc phase (T < Tg) 



(33) 



with the glass temperature Tg ^ C. The corresponding renormahzation group flow diagrams for random 
manifolds and periodic media in a random matrix are shown in Figure El 

In 13 = 1 dimension the roughness exponent for periodic media is C = 1/2 both at zero and non-zero 
temperature ^31^ (compare Appendix Ib)| . 



2.4 Metastability 

The most important feature of pinning is the appearance of metastability. To demonstrate, how metasta- 
bility appears for weak pinning on the Larkin scale, we consider here a D = 1-dimensional example. Let 
us assume that the renormahzation group transformation has been performed n times until the Larkin 
scale Lp — ao5" is reached. Below this scale the perturbation theory is known to be valid and there is no 
problem in deriving an effective Hamiltonian on this length scale. 
The latter can then be written as 

i ^ ^ 

C^"^ and vj(^\i,Ui) are the stiffness constant and the random potential on scale Lp and the Ui are the 
remaining degrees of freedom. By definition of the Larkin scale the elastic and the random part of the 
energy are of the same order of magnitude. In continuing the real space RNG we eliminate now half of 
the degrees of freedom by minimizing the total Hamiltonian, keeping every second (say with i even) 
fixed. This leads to a new effective Hamiltonian 

H^"-^^^= E {^{u,~u,^,r + vj-\^,u.)+SVt-''^ (^,^^^^) \ , (35) 



=2m 



where 



(S\4"+')(i,«) = min [C(")(u - w,+i)2 + Vj;'\i + 1,m,+i)1 • (36) 



The expression in [...] consists of a parabolic potential with the minimum at m^+i = u plus a random 
potential. Since both are of the same order of magnitude, there will be in general several local minima 
(compare with Figure 0)). Let us start in a situation where u = and Wi+i = uj^) is the true minimum 
of this expression. Besides of this minimum in general several other local minima at Ui+i = it(^->x) will 
exist. If we change the external variable u to values different from zero, the minimum Uj^ij may remain for 
small values of u still the global minimum, but eventually for u = another minimum Ui+i = will 

take over the role of the global minimum. At this point the effective potential 5V^~^'^\i^u) is smooth, 
but the derivative jumps by an amount 

2C(")(u(i) - u(„,)) + l42(z + 1, u(^)) - y^;;)(z + l,K(i)). (37) 

As a result, the potential shows a cusp ai u = Uc- Since the magnitude and the positions Uc of the 
jumps of the forces are random, the effective potential acting on the degrees of freedom on scales L > Lp 
is scalloped as shown in Figure ^ Continuing this procedure by eliminating further degrees of freedom 
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Figure 4: Left: The two contributions to the potential Svl^^~^^\i,u) as weh as their derivatives. The 
points of intersection in the lower figure correspond to local minima of the potential. Note that on scales 
L <^ Lp, where the elastic energy is much larger than the random potential, there is only one point of 
intersection and the derivative of the effective potential is continuous. On the contrary, for L ^ Lp more 
and more meta stable states appear. 

Right: The effective potential of Svj^'^^\i,u) as a function of u. The potential exhibits random jumps 
of its derivatives. 

we will obtain more of those cusps on larger and larger length scales. Such a picture was advocated in 
|5f)ll57j : we will demonstrate in the Section El its derivation within the local pinning picture. The forces 
generated by this effective potential V^{i,u) change continuously between two cusps where they jump 
from negative to positive values. The typical distance between two consecutive cusps is of the order I 
since cusps first occur at the Larkin scales. On scales larger than I the forces undergo a random walk 
such that 

{{Vkii^^u) - V^iMf)^ (xl u-u'\ (38) 

for small |m — This structure of the force correlator is in agreement with our RNG analysis. Indeed, 
since —Ruu{u) is the correlator of forces separated by a distance u, the difference 2{R^^{u) — i?Jj„(0)) oc 
— I u I denotes the square of the difference of these forces averaged over the disorder. It has therefore the 
same meaning and the same random walk property as the quantity in equation H38(l |5f>l W>7\ . 

So far we discussed the equilibrium, assuming an adiabatic change of u. If we change however u fast 
enough, the system may not reach equilibrium and remain in the local minimum until this minimum 
disappears completely or it may jump to the new minimum with some delay. Such a situation will be 
considered in the second part of this review when we consider the dynamics of strong pinning. 

3 The close to equilibrium motion of elastic objects under an 
external dc- and ac-drive 

In this Section we want to consider the creep motion of the elastic object in a random environment 
under the influence of a weak external driving force density, / <C /p, where fp = Tp/lLp. (fp is of the 
order of the zero temperature depinning threshold fc discussed in the next Section.) We are now in a 
non-equilibrium situation which requires its own treatment. If / is however small, as wc will assume in 
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this Section, we are sufficiently close to the equilibrium such that we can still use our findings of the 
previous Section (for earlier descriptions of relaxation phenomena in CDWs see e.g. reference |69|). 



3.1 Constant driving force 

We first consider the case of a constant driving force /. All changes in / are assumed to be made 
adiabatically. The equation of motion will be assumed to be overdamped with a bare mobility 7 

- S = - + '?(^' ^ - (x, u) + f + ^(x, t) (39) 

7 ot ou 

where 77(x, t) denotes the thermal noise 

(?/(x, tM^',t'))T = 2-S{^ - K')S{t - t') . (40) 

7 

In the following we will repeatedly consider the elastic object on a variable length scale L which may 
vary from the microscopic cut-off oo to the system size Lq. The coupling between different length scales 
due to the anharmonic random potential in is at least partially incorporated into these considerations 
by (i) the use of the non-mean field exponents Q and % = I? — 2 + 2^ and (ii) the condition, that in order 
to have a moving elastic object on scale L, the system has to be able to move on all length scales below 
L. From the dynamical point of view our present analysis is a type of a mean-field treatment. 

As follows from equations (|18|l . the typical free energy fiuctuations on the scale L are of the order 
F{L) ~ Tp{L / Lp)^ . The energy barriers between different metastable states scale as Eb{L) k, Tp{L / Lp)^ 
with an exponent tp which is in general different from x- In the following we will however assume that 
ip = X- For some systems this can be shown explicitly (see e.g. [TUIE])- In general this is not the case. 
A counter example is a system with an additional isolated potential peak V^5(x — xo)(5(u), which forms 
a barrier which can never be overcome by thermally activated hopping (unless we give up the elastic 
approximation and include topological defects in the structure under consideration). In equilibrium 
statistical mechanics, however, this potential peak does not play a role, since it can easily be avoided by 
the elastic object. This example directs us to the picture of the local pinning considered in the Sec. 
and later on. 

Next we include a small driving force density /. Taking into account that the typical distance between 
different metastable states is of the order w{L) (see Figure ^1, we can write for the expression of the total 
energy barrier 

EB{LJ)^F{L)-fL^wn{L)=Tp(^-^y(^l-(^^^' . (41) 

The second term on the rhs of equation H41|l . ~fL^wu{L), describes the reduction of the barrier due 
to the tilt of the potential by the external force fL^. In rewriting this term we used (|15() and H16|l and 
introduced the force length scale Lj associated with the equilibrium length scale Lp: 

/ f \ 1/(2-0 y 
Lf = Lp[^fj , fp^j^^ClL;\ (42) 

Eb{L, /) is shown in Figure^ it has a maximum aX L ~ Lf = L f{ ^_^_^ )'^/ ^'^'^^ < Lf and vanishes for 
L = Lf. Applying a small driving force / corresponds to testing the system on a large length scale Lf. 
Note, that Eb{L = Lp, /) = Tp[l - (f/fp)] . 
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EB(L,f) 




Figure 5: Energy barrier as a function of the length scale L for a given driving force density /. The 
different lengths and energy scales are explained in the text. 



Assuming for the moment that H41|) . H42|) are valid up to / /p, we see that f ~ fp determines the 
depinning threshold then there is no energy barrier left in the system. Note however that we derived (|41|l 
under the condition / <C /p, the vicinity f ^ fp requires a special treatment which we will consider in 
the following Section. 

For f ,fp the elastic object is restricted in its motion by energy barriers of maximal height 

■lp \ 



EBif)^EBiLfJ)^Tpl^Jj , A^-^, (43) 

which can be only overcome by thermally activated hopping. The creep velocity of the elastic object follows 
from Vcreep ~ ''^{L f) / t{L f) where we use the Arrhenius law for the hopping time r ~ u}~'^e^^'^^^^^ and 
ujp « C7/Lp = jfp/l- This results in a creep velocity 

.(/)«^oce-*(4^)^ (44) 

We omitted the prefactor on the r.h.s. which is beyond the accuracy of the present considerations. This 
formula is valid for T <C Eb{J) and was found first by loffe and Vinokur [72] (see also |73]). Equation 
H44|l was derived from a renormalization group treatment in |74l I75j . In the opposite case T 3> Esif) 
we expect a linear relation between the driving force and the velocity: 

V ~ 7/. (45) 

The border line between the two cases, T « Esif), defines a temperature dependent force fx 

/T = /p(^f] . (46) 
Note that the creep formula is valid only for / <C /t, i-c for f fp and T Tp. 
3.2 Periodically oscillating driving force 

Next we consider the motion under influence of an ac driving force with a finite frequency uj ujp = jfp/l 

fit) ^fo sin (cot). (47) 
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From the Arrhcnius law wc conclude that in driving the system over the period of time ju only barriers 
of maximal height Ei^{T) on a corresponding length scale L^^ given by (compare Figure 1^: 



-E^{T)/T 



1 i.e. E^{T) = Tp 



L,. 



Tin 



(48) 



can be overcome. If E^{T) > Esif), the oscillating force has enough time to trigger jumps over all 
relevant barriers and hence the creep formula (|44|l is still valid (with / replaced by equation (|47|l ). In 
the opposite limit this is not longer the case and hence there is no global motion, v = 0. The relation 
Euj{T) ~ Esif) determines a temperature and frequency dependent cross-over force fu{T) 



1/m 



Tin I 



/t In 



(49) 



which separates the creep region f^j < f fx from the region f < fuj where v = 0. The crossover lines 
fx and fi^ are depicted in Figure El For a discussion of the region f ^ fp see Section^ 

For f < fuj the elastic object as a whole cannot follow the rapidly oscillating external driving field. 
However there is still a local motion of parts of the elastic object corresponding to length scales 



L < = i„ ( ^ In 



1/x 



(50) 



(see FigureEl). We consider to this point in the next Section where we will treat these fluctuating parts as 
two-level systems. As a side remark we note that for similar reasons the specific heat obtains, if measured 

tS- m uj„t 



over a time scale t, an extra factor 



1 



> on the r.h.s. of equation H24|) . 



TA 




v=0 



linear regime 



Figure 6: The cross-over fields /t and fuj{T) as a function of temperature. Note that for w — > fuj 
approaches the T- and /-axis, respectively. For fuj < f < fx the elastic object follows the external field 
in a creep like motion whereas for / < /^^ only segments of elastic object of size L < L^{T) (defined in 
H48|l ') can follow the field and the average velocity of the whole object vanishes. 
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3.3 Dynamic response of the pinned elastic object: Two-level systems 

In this subsection we consider the influence of an external time-dependent field f{t) = f sin cut on a 
pinned elastic object in the region f <^ f^. Thermal motion over energy barriers Eb{L < L^) mediates 
transitions between configurations which have an energy difference AE < T. As we discussed already in 
Section 2.3, the distribution P{AE,L) of AE is smooth and has a width of order F{L) = Tp{L/Lp)^. 
Hence there are only rare pairs of metastable configurations with AE' < T. We therefore model the 
object as an ensemble of noninteracting two-level systems |72l I50j . 

In the following we discuss the dissipation in a two-level system due to an applied ac-field. We begin 
with a discussion of a given two-level systems on length scale L. The separation between the two minima 
is w ~ w{L) and the energy difference is AE. Then the probability that in thermal equilibrium the 
system is in the higher-energy minimum is 

p-AB/T _ 

n„{AE) = ' = [exp (AE/T) + 1] . (51) 



e 



The external field f{t) disturbs the energy difference AE by SEit, L) w f{t)L°WR{L). Therefore the 
system relaxes to the new time-dependent equilibrium configuration n{AE + 5E) w ni^{AE) + 5n{t) . The 
time dependence of 5n is controlled by the relaxation time 

t{L) « cc-p 1 exp {Eb{L)/T), LOp « Cj/LI = jfp/l (52) 

of the two-level system and by the time dependence of 5E{t,L). In a linear approximation the time 
dependence of dn{t) is therefore described by the equation 



d_ 1 



dnp dSE 



The power dissipated in this way by the two-level system of linear size L is given by 

W{L,uj)^-Re(^5n*L''wR{L)^'^ , (54) 

where Re{...)oj denotes the real part of the Fourier transform. With this and the Fourier transform of 
equation H53|l . we get the power absorbed by the given two-level system: 



if / AE ^ "' ~^ 

W(L,uj) cosh 

4T V 2T 



iSEit,L)y ^ ^ , (55) 

1 + W^T^ 



where 6E = fL^WR{L). To get the power dissipated by all two-level systems on scale L...L + SL we have 
to average this expression with P{AE,L), eauation lPT)) . and multiply it with density v{L)dL « ttstt 
(compare with equation (|23|l ) of the two-level systems on scale L. Since the distribution function for the 
AE is smooth and has a width of the order F{L) ^ T, only the fraction T/F{L) of them contributes to 
the average. Hence we obtain for the total power density dissipated by all two-level systems 

/>oo />oo 

Wtotaiiuj) - / dLzy(L) / d{AE)P{AE,L)W{L,uo) 

JLp Jo 

dL 5E^{L) LoML) 



L \LJ F{L) l+tjV2(L) 
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The energy dissipation yV(L,w) is related to the imaginary part of the dynamic susceptibility 

X{L = 27r/ I k |,r,c.) ^ / d{t^t')e^-^'-''\{du^.{t)/df^^{t'))T)B. (57) 



W{L,u)^-u;x"{L,T,u)f\ 



(58) 



The total susceptibility is given by the integral over L with the probability distribution v{L)dL similarly 
to H56|l . The main contribution to the real par x'('^) comes from the length scale L which fulfills the 
condition 9x'(i, w)/9L = 0. This yields (wr)^^ « \xEb{L) /T—\\. For low frequencies and temperatures 
this gives L « L^^ (compare with equation H50|l ) and hence 



x'(T,a;)^ j dLi^iL)P{AE{L),L)x'{L,T,L 



7-2 r 
ZL 
C 



l2/x 



ln(- 



(59) 



Decreasing the frequency leads to an increase of the susceptibility which is a typical experimental trend, 
as well as the logarithmic dependence on the frequency lo. 



4 Critical Depinning at T = 

4.1 Constant driving force, f > fc 

So far we considered the region f fp- Increasing / we expect to reach a critical force density fc (as it 
turns out it is of the order of /p) at which the elastic object is dcpinned. Above the depinning transition 
the elastic object moves even without the help of thermal activation with a finite velocity v which reaches 
the non-critical regime ?; w 7/ at large driving forces /. Qualitatively the v — f diagram is depicted 
in Figure [3 AH changes in / are again to be performed adiabatically. The equation of motion in the 




Figure 7: The velocity of the driven clastic object as a function of the driving force / both at zero and 
non-zero temperatures. 



absence of thermal fluctuations reads 



-^=CW^u + f + g{^,u) = -^-^. (60) 
7 ot ou 
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Here we introduced the pinning force density 

ff(x, u) = -Vr,„(x, u) = — (61) 

ou 

resulting from the random potential. It is also Gaussian distributed with the correlator 

{g{x,u)g{x',u'))R^S{x-x')A{u~u'), A(u) = -i?„„(w) = (62) 

for the bare (unrenormalized) correlators. If the object is completely stiff, i.e. M(x,i) = u{t), ~ this is the 
situation if Lq <C ip, ~, then the average pinning force density {g(yi,u))R vanishes and its fluctuations 
are of the order {Lq^ A{0))^^^ and hence arbitrarily small for a macroscopic object. Thus a rigid object 
would never be pinned. 

Next we shall use perturbation theory for weakly distorted elastic objects To this aim it is 

convenient to go over to a co-moving frame by rewriting u(x, t) = vt + u{'x.,t) with {u)r = and look 
for the lowest non-zero correction to the velocity. Indeed such an approach works well at high velocities 
where the displacements u are small. This gives 

7"^t> - / = (g(x, vt + u(x, t)))R w (gu(x, vt)u{x, t))ii (63) 

where (7„(x, u) ~ dg{x.,u)/du. In the lowest order perturbation theory 

u{x,t) = f d^x' f dt'GQ{x-x',t-t')g{x\vt'), (64) 

J J — OO 

Go(x,i) = 7 / e'kXg-C7k^t 
and with the help of (|62|l . we get from the r.h.s. of equation H63|) 

/>oo 

J~^V - f ^ {AttC)-^^^ dTT~^'^Au{vT/^) (65) 

J(CA2)-i 

The large scale momentum cut-off K ^ a^^ appearing in H64|l was taken into account by adding a factor 
g-fc /A ^ Moreover, we introduced r = i7 as a new variable. We imply here that D > 2 and hence the 
integral (|65|l is convergent at large r. Below we will see that this restriction is unnecessary. Decreasing 
/ and hence v, one finds for v 0+ a non-zero threshold if and only if A„(0+) is non-zero: 

^ _ A„(0+)A^-^ 

(D - 2)C ^ ' 

The result (|66|l looks at a first glance different from fp w CIL~'^ which we found in the previous Section. 
Most importantly, the correlator A(u) = —A(—u) of the random forces has to have a cusp-like singularity 
at the origin since — A„(0-|-) = A„(0— ) has to be positive. This is not what one gets naively from a weak 
random potential Vr(x, u) which shows an analytic behavior of A{u) for small u and hence A„(0-|-) = 0. 

This long standing problem has been overcome by the renormalization group theory of critical depin- 
ning at T = |76l 1771 ITHl 1791 158j . It was shown that the force-force correlator A{u) develops indeed a 
cusp-like singularity on scales L > Lp after the degrees of freedom on scales L < Lp have been integrated 
out. To lowest order in e = 4 — D the renormalization group equation for A{u) in this non-equilibrium 
situation is identical to that for —Ruu(u) following from equation H28f) (calculated under equilibrium 
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conditions) by differentiating R{u) twice with respect to u. Note that this simple relation breaks down 
to order [SHI- Hence A,j(u 0) 7^ 0. The force correlator on these scales becomes scale dependent 

A(^.) ^ A(u; L) oc {Cl/Li)'L-^+''+^iA* {ti{L/ L.p)-~^ /l) , (67) 

where the function A*(y) has a cusp for small y. Similarly, the mobility 7 is replaced by an effective scale 
dependent expression 

/ r \ 2-i 

7 -> 7(i) « 7 (^;^ j • (68) 

C and z are two new non- equilibrium critical exponents which can be calculated by an expansion in 
e = 4 — _D. Note that these exponents are in general different from the equilibrium exponents introduced 
in the earlier Sections. If we replace by Lp and A„(0) by A(w, Lp) in equation (|66(l . and put C = 
(since the corresponding integral is dominated by small scales L w Lp), we indeed arrive &t fc ~ fp- A 
detailed calculation gives |Z3 EH! 

fc = — ^A;(0+)^. (69) 

Qualitatively, the cusp singularity can be understood as follows: In order to obtain a non-zero depinning 
threshold the average value of lim,„_,o(5(x, vt + {((x, t)))^ in equation H63|) has to be negative. In other 
words, the elastic object "sees" in a pinned configuration more increasing than decreasing potential hills, 
even at u — > 0! ^ The next step in (|63|1 is to expand g(x, vt + u) with respect to u. This leads to a 
product {gui^, i^^).9(x, vt')) ji with t > t' , which, in the limit v ^ (hence vt uq, vt' — > mq — e„) has to 
be negative as well: 

(,9„(x, uo)5(x,Mo - eo)>H ~ (.g„(x,Mo)[5(x, uo) - et.,9«(x, uo)])_r < (70) 

For a typical potential dominating the correlations of the random forces two cases are possible: 

(i) If 5«(x, uq) < 0, i.e. the force is locally decreasing with increasing m, then the force at uq — e„ has 
to be positive, i.e. accelerating. This is the situation shortly before one reaches a potential minimum. 
For 5 < one concludes from the r.h.s. of (|7n|) |5m(mo)| > \9{uo)\/£v, i-C- in the hmit v ^ ey ^ the 
curvature of the potential of pieces with g < becomes arbitrarily large and correspondingly these pieces 
of the potential shrink to zero. 

(ii) If on the contrary (7„(x, uq) > 0, i.e. the force is locally increasing with increasing u, then the 
force at uq — has to be negative, i.e. retarding. This is the situation shortly before one reaches a 
potential maximum. For g > one concludes from the r.h.s of H70() gu{uo) > g{uo)/ey, the curvature in 
these pieces of the potential diverges for ?; ~ e„ ^ and hence these pieces of the potential disappear as 
well. 

Therefore for v ^ the effective potential consists mainly of pieces where g and gu have different 
signs. Regions with the same sign of g and p„ disappear gradually from the effective potential emerging 
on scales L > Lp (or at least give only a small contribution to it). If we assume that only pieces with 
ggu < indeed remain in the effective potential and assume that these pieces extend up to ggu = 0, 
then a continuous and piecewise differentiable potential can be constructed from alternating segments 
where g and —gu are both negative or positive, respectively. Such potential pieces are given by V+(u) for 

^As mentioned already, within the collective pinning regime this is only possible on scales larger than the Lar kin-length, 
since for L < Lp the elastic object is essentially undistorted. More generally, the cusp appears together with metastable 
states; it follows from the possibility to switch between ascending and descending branches (in terminology of the Section 
17.11 preferably selecting the lowest one (see the derivation of the force correlator in the Section Il2l . 
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< w < w+ and V- (u) for u_ < u < which wc can model as e.g. (compare Figure |SJ: 

V+iu) = Uu (l ^ V-(^) - /-— fl - I—) ■ (71) 

The forces f± and positions u± may change from segment to segment but have to guarantee the con- 
tinuity of the potential. According to what we said above, the average value of ~ 0{fp) has to be 
identified with the depinning threshold. Such a shark- fin potential is schematically drawn in Figure 
Its appearance in a driven situation is rather obvious since the elastic object will pass the regions where 
gug > very quickly. Together with the property A(w) = A(— u) this explains the physical origin of 
the cusp. (See also the above footnote; the precise meaning of the "shark- fin" potentials is given by 
"termination points" of metastable branches of the Section l7. II and corresponding figures.) 
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Veff(u) 
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/ \ 















Figure 8: (a) The effective potential as it emerges on scale L > Lp for an elastic object moving with 
a positive velocity. The shape results from the requirement that A^(O-I-) is negative and that at the 
boundaries of each segment gu{u)g{u) ~ 0. (b) The effective potential resulting only from the condition 
that gu{u)g{u) is negative everywhere. We consider the case (b) to be the generic one. Note, that 
orientation of the potential is reversed for a motion with a negative velocity. 

For a negative driving force a corresponding discussion leads to a potential in which g and gu have the 
same signs, regions with the different signs of p and gu disappear. This corresponds to a potential which 
differs from that drawn in Figure |H1 by changing the direction of the u- axis. If we do not impose the 
additional condition that ggu = at the boundaries of the intervals where gu > or gu < 0, respectively, 
then the potential may be discontinuous as depicted in Figure 

Another problem following from perturbation theory is the fact that for D < 2 the r.h.s. of equation 
l)65|l could diverges. This also has been overcome by the renormalization group theory of the critical 
depinning |76l 1771 17?^ I58| : = 2 is not longer the lower critical dimension due to the appearance of 
non-classical critical exponents in the renormalized perturbation theory. 

Close to the depinning transition the velocity - which can be considered as an order parameter of the 
transition - vanishes as a power law 

, />/c. (72) 
Here we introduced the characteristic velocity scale Vp = jfp- Approaching the depinning transition 



f f-fc \ 
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there is a diverging correlation length 



f-fc 



fc 



(73) 



The appearance of a diverging correlation length on both sides of the depinning transition has to be 
expected for the following reason: if we approach from values / < /c by increasing / adiabatically, 
larger and larger avalanches of local motion of the elastic object will occur until we reach a critical state 
at f = fc- A further increase of / will then lead to a macroscopic motion of the elastic object. At f = fc 
a local perturbation will hence trigger a global response, corresponding to an infinite correlation length. 
With increasing velocity spatial fluctuations in the local velocity will be reduced and the correlation 
length shrinks again. We will come back to this point in the following. 

On length scales ^ L <C ^ the non-linearities of the pinning potential dominate and distortions of 
the elastic object obey dynamical scaling with non trivial exponents C and z: 

((^(x,t) - u(x',0)'>f = I (^^)'* ((^^)7-p(^- . (74) 

Here C is the non-equilibrium roughness exponent and ujp ~ Vp/l. It turns out that 1 < z < 2, i.e. the 
dynamics close to the depinning transition is super-diffusive, reflecting the rapid motion of the object 
after the maximum of the shark fin potential has been overcome. The scaling function $(y) behaves as 
y-C/^ for y Q and approaches a constant for y oo. The critical exponents satisfy the new scaling 
relations |77] 

' ' ' (75) 



2-C S-C D+C 

These exponents were calculated first to order e = (4 — _D) in [77j and recently to order [SS]. For charge 
density waves C = |7S](i.e. the roughness increases logarithmically with L) and z = 2 — | — ^ [7^155] . 
whereas for domain walls ( ^ ^{1 + 0.14331e) and f = 2 - f - 0.04321e2 

In the opposite regime L ^ ^ the problem is essentially linear and u can be replaced by vt in 
the argument of g{'x,u). This can be seen qualitatively as follows: On the time scale t the elastic 
object advances on average by an amount vt. Randomly distributed pinning centers will lead to a local 
distortion which, according to (|7^ . spreads over a region L(t) « Lp (tOpt)'^^^. The local retardation 
or advancement of the object due to the fluctuation in the density of the pinning centers scales as 
u(t) « l{L{t)/Lp)^ « l{ijjpt)^/^. Since ^ < z, u{t) grows more slowly than vt. Thus on time scales 

t > ty = LUp^ {vp/vY^^''~'^^ and length scales L > S, = L(ty) the non-linearities in the argument of 
Vr,u{^, vt + u) can be neglected and the Hnearized theory applies. In this case m(x, t) is replaced by vt 
in the argument of the random forces. Random forces act then as a thermal noise with the temperature 

^v-K 

So far we considered the elastic theory of critical depinning. If we include topological defects in the 
theory, the transition may become hysteretic, as was shown in [SO]. We will come back to the influence 
of topological defects on pinning phenomena in the second part of this article. 

As a side remark we mention here that an alternative characterization of the depinning transition can 
be reached if we pin the elastic object at the boundary of the system by an infinitely strong surface barrier 
such that the displacement at certain surfaces vanish. In a charge density wave this can be reached by 
applying an external electric field but preventing a current fiow by the absence of external leads. In a non- 
random elastic system an external force / then generates a parabolic displacement profile as a solution 
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of (|60|l . The situation is different in a system with random pinning forces: as long as / < /c the elastic 
object cannot move and the pinning by surface barriers does not matter. For / > /c on the other hand 
the surface barriers prevent the elastic object from moving and a parabolic profile will emerge. A detailed 
investigation shows that this is indeed the case |SI]. Using the decomposition ii(x) ~ uo{x) + u(x) with 
the Ansatz 

1 ^ 

uo{:x.) = --^Ci{xi - Xifi)"^ + uo (76) 

4=1 

and {u(x))ii = 0, one can determine fc from the vanishing of C = X^i^i- Which of the curvatures Ci 
are non-zero depends on the specific pinning conditions on the surface. Note, that this decomposition 
is similar to the description of the dynamics of the depinning transition in a co-moving frame. From 
equation H6()() we obtain 

CC-/ = (.g(x,wo(x) + H(x)))«(54x,uo(x))u(x)) , (77) 

which replaces equation Ht)3|) of the case of a moving elastic object. In lowest order of perturbation theory 
we obtain then (similar to the derivation of equation H65|) ') its RG-counterpart [51] 

C = ^^. (78) 

It is to be expected that this relation is true to all orders in e = 4 — D. Indeed, on the scale of the 
correlation length ^ the height a^^ of the parabola is expected to scale like the roughness 1{£,/Lp)'' . Thus 

and with (|73() . \T&\ and the scaling relation equation (|75() v = we get indeed the expected result. 
Thus we may also characterize the depinning transition by the vanishing of the parabolic displacement 
profile. If one decreases the forces again the curvature shows a pronounced rhombic hysteresis profile 
[5T] . The problems of inhomogeneous profiles are related to contemporary space resolved studies of sliding 
CDWs, see the Section El for discussion and references. 



4.2 The depinning transition at finite temperatures 

At / < /c and T = the velocity is zero, but one has to expect that as soon as thermal fluctuations 
are switched on, the velocity will become finite. Scaling theory predicts in this case an Ansatz |H21 155] 
(generahzing ifT^ ) 

v{f,T)^TP'^^{L0^ (80) 

with $(a;) const, for a; ^ and $(x-) ~ for > 1, such that v{fc,T) ~ T'^/^. Here r > is a 
new exponent which still has to be determined. 

This prediction seems to be in contradiction with simple scaling considerations applied directly to the 
equation of motion l|39(l . Indeed, after renormalization of this equation at T = up to length scale ^ and 
time scale uj~^{£^/Lp)^ all terms in this equation scale as (/ — fc)- If we consider now the thermal noise 
as a small perturbation (at the fixed point describing the depinning transition) then its contribution to 

^ / t\^^^ / f-f 

the equation of motion is of the order [f — fc) j ( j j , where x = D — 2 + 2(^>0 and hence 
r = — i'x < 0. Thus thermal fluctuations seem to be irrelevant at this transition. 
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However, this is not true. The previous argument considers the influence of thermal fluctuations on 
length scales of the order L w ^. The relevant thermal fluctuations which depin the elastic object act 
however on much smaller scales of the order of Lp ^ as we will see now, following an earlier argument 
by A. Middleton At the critical point / = fc essentially only barriers on the scale L k, Lp are left as 
we saw in the previous Sections. It is therefore sufficient to consider only this length scale. To this aim 
we coarse-grain the system into regions of linear size Lp and denote the corresponding Larkin-domain 
by the subscript i = 1, TV = {L/Lp)^ (compare with Figure EJ. 




Figure 9: The decomposition of the system into different Larkin domains (left), the effective potential 
for the coordinate Ui in the case of a smooth (middle) and a ratcheted potential (right) 

Each domain is then essentially characterized by a single local coordinate Ui related to w(x) by 
Ui ~ J^d^ xu{x.), where j^d^x denote the integration over the «-th domain. Treating the interaction 
between different domains within a mean-field approximation one can write (in spirit of the Kim- Anderson 
approach) a local energy expression for the domain i 

n^^V,{u,)-fu,^n,{f-f,)u, (81) 

Vi{ui) denotes the effective potential for the coordinate Ui at the depinning threshold f = fc- The term 
—Ki{f — fc)ui, Ki > describes the mean-field type coupling to neighboring domains. The full energy is 
then given hy H = Y^fL^Ki. 

Two types of effective potentials were considered j83) : (i) a smooth potential Vs{u) = /pu(l — 
and a ratcheted kick potential with Vrkiu) = /pu(l + (u < 0, compare with Figure|5J). uq denotes the 

position of the potential minimum. In general the values of uq (> 0) and fp (> fc > 0) will be different in 
different domains. For simplicity we assume in the following that uq ~ ILp is the same for all domains. 

At the threshold / = fc there is a metastable state corresponding to the left minimum of the potential. 
The most unstable domain is then characterized by a value of fp which slightly exceeds fc and « 1. If 
we increase now / from f = fc 'to f = fc + S f , then the minima in all domains with fp — fc < + Ki)6f 
will disappear. These unstable domains will trigger transitions in neighboring domains which destabilize 
further domains and so on until the whole object is depinned. The height of the energy barrier for the 

most unstable domains in the region Sf < is given by Eb ^ u^fc ( (1 + k) y 1 for the smooth 

and Eb ^ uofc (^(1 + k) ^"^^ ^ fo'^ ^^le ratcheted kick potential, respectively. These barriers become 
irrelevant at temperatures T > Eb- An increase of T from T = has the same effect as increasing / 
by a value ~ T'^/'^ from which we conclude r = 3/2 for the smooth potential. For the ratcheted kick 
potential the increase of T has the same effect as an increase of / by a value ^ T^/^ and hence t = 2. 
The different exponents reflect the different non-linearities of the effective potential. Our analysis of the 
effective potential depicted in Figure |Hla (with the analytical form given by equation l|7T|l ) is different 



28 



from |83| . but the exponent t = 3/2 is identical with that of the smooth potential, since the barriers in 
both cases have the same dependence on (fc — f), as one can check easily. However this remark is only 
true for a potential of the type depicted in Figure |Hla, for a potentials of the type depicted in Figure 
|Slb the exponent r is probably non-universal. The exponent t found from simulations for domain walls 
in random field systems jM] varies indeed in the range 1.5 < t < 2. Naturally the exponent t = 3/2 
appears in treatment of termination points fSections I7.1IC|I which specify the meaning of the suggested 
above instabilities. 



4.3 Depinning due to an ac-field 

In this subsection we want to discuss the effect of an ac-field of a finite frequency lu on the pinning of 
the elastic object at zero temperature. The equation of motion is still given by H6()() with / — > f{t) = 
fsm{ujt). A finite frequency uj of the driving force acts as an infrared cutoff for the propagation of 




Figure 10: Schematic frequency-field diagram for the depinning in an ac external field (with / > fc): 
For < cij <C the depinning transition is smeared but traces of the to = transition are seen in the 
frequency dependency of the velocity at / = fc- This feature disappears for lu ujp. 

perturbations, resulting from the local action of pinning centers on the object. As follows from H6U|) with 
the renormalization H()8|) these perturbations can propagate during one cycle of the external force up to 
the (renormalized) diffusion length 

= Lp(7C/c.L^)i/^' = Lp(c.p/c.)l/^ = ^ . (82) 

If Lf^ < Lp, i.e. to > Up, then there is no renormalization and z has to be replaced by 2. During one 
cycle of the ac-drive, perturbations resulting from local pinning centers affect the configuration of the 
elastic object only up to scale L^, such that the resulting curvature force CIL^^ is always larger than 
the pinning force - there is no pinning anymore. 

In the opposite case L„ > Lp, i.e. uj < ujp, the pinning forces can compensate the curvature forces 
at length scales larger than Lp. As a result of the adaption of the elastic object to the disorder, pinning 
forces arc renormalized. This renormalization is truncated at L^^. Contrary to the adiabatic limit uj —f 0, 
there is no sharp depinning transition if w > 0. Indeed, a necessary condition for the existence of a sharp 
transition in the adiabatic case was the requirement that the fluctuations of the depinning threshold in 
a correlated volume of linear size ^, 5fc ~ fc{Lp/S)''-^^^^/'^, are smaller than (/ — fc), i.e., [D -\- > 2 
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(compare equation H75|l ') )77| . For uj > the correlated volume has a maximal size X^; and hence the 
fluctuations Sfc are given by 



Sfc 
fc 



^ X {D+0/2 / ^ X (15+0/(25) 



(83) 



Thus, different parts of the elastic object see different depinning thresholds the depinning transition 



0.4 
0.2 



^ 



-0.2 
-0.4 



fo=0.8, 00=0.3 




fg=f^ , a)=0.3 




03=0.0 

-fc 


.i.^----'T--y'''''.-^-^^ ^^^^^^^ 






[ f~c fc 



-0.6 -0.4 -0.2 




f(t) 



0.2 0.4 0.1 



Figure 11: Velocity hysteresis of a D = 1 dimensional interface in a random environment. 

is smeared. Sfc has to be considered as a lower bound for this smearing. A full understanding of the 
velocity hysteresis requires the consideration of the coupling between the different L(^-segments of the 
elastic object. Approaching the depinning transition from sufHciently large fields, f{t) ^ fc (and w ^ ujp), 
one first observes the critical behavior of the adiabatic case as long as ^ <^ L^. The equality ^ ~ 
defines a field fco signaling a cross-over to an inner critical region where singularities are truncated by 

Note that fco — fc ~ fd^/uipy^'^''^^ > Sfc (cf. FigureEB- I* is then obvious to make the following 
scaling Ansatz for the mean interface velocity (/o > fc, Vp = ujpl) 

Here the subscript ± refers to the cases of / > (<)0, respectively, and (p±[x ^ oo] ^ . For f{t)~fc ^ fc 
the classical exponent (3 — 1 applies. For \x\ <C 1, <y9± approaches a constant c±. Function ip- changes 
sign at a critical value fd^j) ~ /c(l — C-{w /iOpY^''^'^^). The velocity shows a typical double hysteresis 
fFigure [TT|l . Qualitatively, the hysteresis loop can be understood to result from the motion in the ratchet 
like potential. Figure |H1 The reader is referred for details can to reference There is interesting 

related work on the influence of an alternating current on pinned vortex lattices by Kohandel and Kardar 

ESIIHZI. 



m 

fc 



(84) 



5 Macroscopic perturbations and external constraints. 

Now we shall discuss briefly the influence of external forces or constraints on the statistical properties 
of our model equation (Q) . The topic includes the important case of topologically non-trivial distortions 
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which can be enforced by appropriate boundary conditions or appUed external forces. We will assume that 
these constraints are kept constant or are changed only adiabatically such that we can apply equilibrium 
statistical mechanics. We will assume that a field CA(x) couples in the Hamiltonian linearly to Vm, i.e. 
there is an extra piece STi.A h^ the Hamiltonian 

SHa^ J d^xCA{x)Vu. (85) 

Examples for A are given e.g. by 

(i) a constant external force / coupling to u for which CA = fx/D or 

(ii) by a field which enforces a dislocation line into the system (if we consider a periodic systems - see 
the following Sections). In the latter case A obeys the relation <f^ Adx = —b = 2TrnQ where the curve C 
encloses a dislocation line and n is integer. 

In general, extra pieces of the Hamiltonian of the form (|85|l lead to an instability: the system feels 
a constant driving force or a torque. In order to prevent an unlimited response we have to assume the 
existence of additional surface forces which keep the system in equilibrium. This will be done in this 
Section (for more details see 

It is convenient to go over to the new field u by 

u(x) = w(x) + / A(y)dy. (86) 
Jo 

In the case (ii) u(x) may depend on the path along which the integration is performed. Different pathes 
may lead to changes of u(x) by mQ where m is integer. Since the random potential is periodic in u with 
periodicity 2Tr/Q such an ambiguity is however irrelevant. After this transformation the Hamiltonian is 
rewritten as 



H + SnA = I d^x I ic(Vu)2 + Vr (^x, u(x) - ^ A(y)dy^ - ^CA^j 



(87) 



In both cases (i) and (ii) the new Hamiltonian H87|l has the same statistical properties as the original 
one, equation since Vr(x, u(x)) is a random function of both arguments. This can most easily seen 
by using the replica method, in which the disorder averaged free enthalpie 



(G{A})„ = - /lnrre-(«+*'«^)/^\ = - lim - 



Tre^— - 1 



follows from the replica Hamiltonian 

Hn = fd^'xY,^ {{"^UafSa.b - ^R{Ua ' "h) - nA^ . (89) 
a,b=l ^ 

Apparently, the replica Hamiltonian is unchanged, apart from the additional term —n J d^x^A^ . It is 
worth to mention that this is true only if the random potential Vr(x,u) is strictly uneorrelated in x. If 
the correlations are given by a smeared out i5-function of width oq, R{uai^) ^Ufc(x)) in equation l|89|) has 
to be replaced by / d^x'R[ua{x.) - M6(x') - A{y)dy]Sa„ (x - x'). 
The disorder averaged free energy follows then as 

(n(V«)})^ = (G{A})^ - J d^xCA{Vu) = (^^{0}) + j d^x^{Vu)lR, (90) 
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where 

{Vu)t,r = C-'^SG{A}/5A = -A. (91) 

If A represents a dislocation then also the mean displacement {u)t,r shows a dislocation structure. 
Correlation functions of m(x) in the presence of external forces can now easily be calculated by using the 
decomposition equation (|86|1 . since u is not affected by the presence of A. 

The glassy phases discussed previously have been found under the assumption that topological defects 
have been excluded. Wc may now consider the stability of these phases with respect to topological defects. 
In particular, we will briefly consider the stability of the Bragg glass in charge density waves phase with 
respect to dislocations. 

Adding a dislocation increases the disorder averaged free energy according to eauation H9()(l by / d^x^A^ = 

l^i^^^ln^, Lq denotes the size of the system. This expression is the energy of the dislocation line 
(D=3) in a pure system and hence dislocations seem to be always disfavored. However the dislocation 
may take advantage from fluctuations in the disorder distribution and choose a position where its energy 
is lowered with respect to the average value. To this aim one has to consider the sample to sample 
fluctuations of the free energy (_F^{(Vu)})_r — {F{{'^'^)})'r- This is a difficult problem and only prelim- 
inary results exists, which support the existence of a quasi-long range ordered phase in d = 3 dimensions 
|62[ 1631 IMl I65| provided the disorder is sufficiently weak. In c? = 2 dimensions, where the disorder in the 
forward scattering term growth under renormalization as InL, dislocation always appear |88|. 

6 Plastic deformations and topological defects. 

The starting point of the collective pinning picture considered so far was that the displacements u grow 
unlimitedly at large distances. At the same time, local deformations (i.e. strains - gradients Vm) were 
assumed to remain small thus allowing for the universal elastic media description, expansions of bare 
energies in terms of Vu or forces in terms of m, etc. Nevertheless the elasticity can be broken already 
at the local level in which case we refer to plastic deformations [3^1 . The effect is not only quantitative 
which would simply affect basic parameters. It happens that plastic deformations related to impurities 
can originate metastable states which is the principle ingredient of the pinning picture. These plastic 
metastable states create a set of pinning effects of their own nature, but also they clarify, or even challenge 
sometimes, the complex picture of the collective pinning (see more in the Sec. I12II . 

In general, plastic deformations invoke displacements which are not small at a microscopic scale, e.g. 
the domain wall width or the crystal periodicity. The deformations may be topologically trivial like large 
curvatures of domain walls, or as vacancies and interstitials in Wigner or vortex crystals; then plasticity 
comes from the strong pinning potential itself. The plastic deformations can be topologically nontrivial, 
and these are locally stable even without impurities. Among our cases they appear only in periodic 
systems because of their ground state degeneracy x x -|- b, where b is any of d primitive periods of 
the sliding crystal. Here the topological defects acquire forms of dislocation lines, or dislocation loops 
and their particular limits of solitons in quasi Id systems, (concerning dislocations, see jSH] for a general 
review, [361 19(J| for the theory, |^|^ for the special case of CDWs and for helpful illustrations.). 

The basic object is the dislocation line crossing the whole sample (reduced to the dislocation point in 
d = 2 systems), see the Figure^] It is a kind of a vortex of displacements u (which are now the D ^ d 
dimensional vectors u =J> u) such that going around the dislocation line u acquires a finite increment 
i5u called the Burgers vector. In principle, it can be any allowed translation of the regular lattice, but 
only dislocation lines with minimal values of 6u, coinciding with one of the primitive translations (5u = b 
are stable which we shall imply below. In our perspective, finite displacements only along the sliding 
direction x are important, so in our studies we shall assume that b = (6, 0, 0) is chosen while dislocation 
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Figure 12: Topological defects in a CDW. The solid lines describe the maxima of the charge density. The 
dashed lines represent chains of the host crystal. From left to right: dislocations of opposite signs and 
their pairs of opposite polarities. Embracing only one chain of atoms, the pairs become a vacancy and 
an interstitial or ^2tt solitons. Bypassing each of these defects, the phase changes by 2tt thus leaving the 
lattice far from the defect unperturbed. 



lines and dislocation loops lie in the perpendicular plane, these are the so called edge dislocations. In 
d — 3, the dislocation line cither must cross the whole sample or it must be closed to the dislocation loop 
(two dislocation points D,D* in d=2). All paths across the dislocation loop acquire the displacement 
b in compare to outer paths. The smallest dislocation loop embraces just one line of atoms with one 
unit cell missed or acquired along this selected line. This limit is the ±2tt soliton in quasi Id systems 
(for short reviews on theory see [301122 and also [HI], for experimental aspects see [HH])- In isotropic 
crystals (Wigner crystals, vortex lattices) the elementary dislocation loop is the symmetry broken state 
of the vacancy or the addatom. Going along two paths parallel to b, one above and another below the 
dislocation line, the difference b of lattice displacements will be accumulated. Then the dislocation line 
can be viewed as the leading edge for an additional atomic plane being introduced to (or withdrawn 
from) the crystal (from the side boundary, or from another dislocation line - the counterpart D* with the 
opposite circulation b* = — b). 

Here we already arrive at the first general significance of dislocation lines for sliding crystals: their 
necessity to bring in or modify the sliding regime providing the so called phase slip processes^. Within 
the CDW language, the phase slippage is required at junctions for the conversion from free to condensed 
carriers |9f)ll97ll^ . When the CDW is depinned between current contacts, CDW wave fronts are created 
near one electrode and destroyed near the other, leading to CDW compression at one end and to its 
stretching at the other end. In a purely ID channel, the order parameter can be driven to zero at once 
|96j which allows for the macroscopic phase slip. For samples of finite cross-section, phase slippages 
develop as dislocation lines proliferate across the sample, each dislocation line allowing the CDW to 
progress by one wavelength |^. Proliferation of dislocation lines or expansion of dislocation loops is 
called the climb. As for any motion not parallel to b, the climb is not conservative with respect to 
the number of atoms (the charge in electronic crystals). As such, it is ultimately related to the current 
conversion requiring for the phase slips. 

Effects of phase slips and the current conversion are closely related to macroscopic strains of the 
sliding and/or pinned state. Recent years have brought a new understanding of the fact that the sliding 
state is also essentially inhomogeneous |1()21 11031 ll()4| . A freedom for deformations is demonstrated by 
the dilemma associated with the choice of a solution for the generic equation describing the sliding motion 
and deformation of the CDW phase: j'^dtf + Cd^tp ~ f (see the Section 4). Taken alone, this equation 
is satisfied by any solution of the type: (p — f{cjt+ (1 — c)a;^/2C) with an arbitrary value of the partition 
coefficient c. Then, at first sight, the response to the driving force / is optionally distributed between the 

^ Phase slippage is a common phenomenon in condensed matter systems with complex order parameters. It has been 
intensively studied in narrow superconducting channels 1991 llOOl . in superfluid helium IIOII and in quasi one-dimensional 
CDW systems. Phase slips have been incorporated to the picture of the collective pinning only recently 1681 



33 



viscous ~ t and the elastic ^ reactions, leaving the collective current undetermined. It is specifically 
the equilibrium with respect to phase slips which selects the solution c = 1 leaving only the viscous 
nondeformed regime ~ t in the bulk [1041 \WS\ . Phase gradients originating from the external force 
cannot grow indefinitely because the associated strain is released through 27r-phase jumps repeating in 
time, with a rate dependent on the magnitude of the remaining strain 
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Figure 13: Generation of the perpendicular flow of dislocations by sliding through narrows. The encircled 
up-arrow indicates the proliferation of a dislocation line. The up-right double arrow indicates the material 
flow; its conversion provides the climb of dislocations which results in the velocity enhancement. 



There is a particular case of the conversion with the help of dislocations which brings us closer to the 
problem of strong pinning. Namely, consider the sliding along a step-shaped host sample (which is quite 
an important issue in reality) as shown on the Fig. 1131 Coming to the threshold, the sliding in cut layers 
is terminated, so that in the narrows the phase velocity v = —ip must increase f — > u + Av. To keep the 
crystal connectivity, new periods must be introduced with the phase slip rate Av/2tt which are provided 
by the flow of dislocation lines in the cross-section. Understanding this macroscopic example leads us to 
expectations for the role of plasticity and topological defects in the pinning problem, even beyond the 
strong pinning limit. 

Let us reduce the size and the sharpness of the obstacle (remind soft macroscopic defects studied by 
the space resolved X-ray diffraction |1(J7| '). That will be then a local region of the enhanced pinning force. 
It can be cither mesoscopic, originated by rare fluctuation of the collective pinning strength^ (e.g. the 
concentration of impurities), or microscopic: a single strong pinning center. Our only requirement is that 
the local pinning enhancement is strong enough to reduce the mean sliding velocity that is to provide - 
at least from time to time - the retardation by the whole lattice period. Then the retarded zone must 
be surrounded by dislocation loops to provide the matching with the rest of the crystal. 




Figure 14: Motion of a sliding atomic array through the strong attractive pinning center, as viewed 
from the co-moving frame. The upper row: the straight arrow indicates the relative displacement of the 
impurity together with the trapped atom. The lower row: having been displaced by more than half of 
the period, the trapped atom is released to its, distant now, position (the left arc arrow), while the more 
close now next atom is trapped instead (the right arc arrow). 

^Thc phase slip rate is given by the space- time vorticity / = {dtdxf — dxdtip) /2tt. 

''interpretation of pinning in terms of large scale fluctuations of the pinning force was suggested in 11061 . where also the 
phase slips have been discussed 
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The microscopic case is illustrated on the Figure ^1 for an attractive impurity moving across the array 
of "atoms" of e.g. the Wigner or the vortex crystal. Usually we shall assume the co-moving frame where 
the pinning center moves through the asymptotically immobile crystal. 

There are three apparent regimes. A weak attractor will only perturb atoms which will smoothly 
return to their equilibrium position. A stronger attractor will draw the atom for more than a half of 
the period, then it becomes apparently favorable to release the overdrawn left atom and catch instead 
the next atom at the right which is now closer, as it will release the energy of deformations. Finally, 
let the attractor be strong enough to draw the initial atom over the whole period to the next regular 
atomic position. Then the crystal comes again to the local equilibrium but in expense of creating zones of 
dilatation and compression behind and ahead of the attractor. Being integers of the atomic period, these 
deformations correspond just to dislocation loops of the minimal size, the solitons of the CDW language, 
embracing the path of the attractor motion, see the Figure IT^ Their energy 2Es will be payed for the 
pinning preserved over one period. 
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Figure 15: A pair of dislocation loops generated by a strong pinning center after the nearly complete 
period of sliding. The cross-section (the figure's plane) corresponds to the quaternion of dislocations in 
the r.h.s. of the Figure El The phase deficit at the impurity point Xi is determined by the steric angle 
of the loop. 

An opposite case to the atomic lattice is the CDW which density is smoothly distributed over the 
whole period, see the Figure El Now the development of elastic and plastic deformations in the course 
of the impurity motion is illustrated on the Figure [TCI with the same consequences for metastablc states. 

=0+2solitons 

271 
71 
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Figure 16: Evolution of the phase profile Lp{x) (for the chain passing through the impurity) in the course of 
the relative motion. Starting from the equilibrium position when 1^9(2;) = 0, the profile evolves gradually 
through the shapes 1,2,3,0 finally developing the bisolitonic shape. These configurations correspond to 
the retarded branch i?+ which becomes metastablc after (f{Q) crosses tt. Since then, the advanced profiles 
l',2',3',0 of the branch E- are less deformed and hence cost a smaller energy W. If the relaxation 
E+ E- does not happen, the new circle starts with the profile = -|- 27r corresponding to the 
infinitely divergent pair of solitons. For a weak impurity the level tt is never reached and only a smooth 
reversible evolution is allowed following shapes 0, 1, 0, 1', 0. 

In summary, we can generally anticipate several regimes which existence can be verified for particular 
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models considered below. We shall refer to some parameter V characterizing the pinning center strength 
with respect to certain thresholds Vi, V2 for plastic deformations. 

1. V <V\. The pinning center is very week, local deformations do not grow to the plastic threshold 
but they smoothly return to the original unperturbed state after the whole period is passed in the course 
of sliding. 

2. Vi < V < ¥2- The pinning center is strong enough to provide a retardation for more than the 
half of the sliding period (|A(p| > tt in the CDW language). Since then the branch becomes metastable: 
it is favorite to switch the deformation from the overdeformed retarded configuration — 27r < Aip < 
— TT, \Aip\ > TT to the weaker deformed advanced one < A(p < n, \Aip\ < tt which saves the energy of 
elastic deformations around the pinning center, see the Fig. A coexistence of stable Ei and metastable 
E2 branches - the absolute and the local minima of the energy - implies the existence of the third branch 
£"3, the energy maximum, which is the barrier separating the two minima. Since we postulated here that 
the metastable branch cannot be maintained over the whole period, then there must be a termination 
point (the "end" point) 6^ where the metastable and the barrier branches merge together to disappear 
at larger retardation, see the Figure [T7I 

3. V2 < V . An even stronger pinning center can sustain a retardation by the whole period or more, 
then there is no termination point and all branches always coexist. After the whole period passes, the 
pinning center is again at the equilibrium with the surrounding crystal, but it is the equilibrium modulo 
27r, the difference of one period is accumulated between the pinning center vicinity and the crystal at the 
infinity. And now this is just the job of dislocation loops to compensate for this difference. Their diverging 
pair forms a cylinder (containing the pinning center) where this difference is just provided. In the quasi- 
Id picture, see the Figure El this process is easily visualized as a creation of the soliton-antisoliton pair 
(the bisoliton) which opens the 2tt retarded allowed ground state in between. 

i 

metastable E2 

stable El 

J 

Be 

Figure 17: Energy branches for a restrictedly bistable impurity. The uppermost thin line shows the 
barrier branch E^. Thick lines show the locally stable branches E±, also classified as E2 > Ei. The 
difference AE ^ E2 — Ei gives the dissipated energy. The difference U = E3 — E2 gives the activation 
energy for a decay of the metastable state i?2. 

Notice finally that the " mesoscopic" case of a cluster of impurities will have more degrees of freedom 
which can originate a large number of close metastable states, thus merging gradually with the collective 
pinning picture. 

7 Local metastable states. 

All the following content will exploit the efficient language of phases for Charge/Spin Density Waves 
(CDW/SDW), which is widely used since ^U]. The CDWs are characterized by the sinusoidal density 
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profile ^ cos(Qx+<y5), sec ((TJ, and have elastic properties of uniaxial crystals, see the Table 1. The order 
parameter can be taken as ~ exp[i<p] so that dislocations are easily viewed as usual vortices. For a periodic 
sliding media in general, the natural choice for the microscopic length scale is the unperturbed lattice 
period b along the sliding direction. It corresponds to the CDW convention to use the phase for the 
description of the displacements u ~(p: ip = ~2Tru/b = — Qu. The velocity becomes the phase velocity 
for which we shall use the same notation v = du/dt ^ —dip/dt. This phase velocity is accessed directly 
in experiments by measuring the so called Narrow Band Noise (NBN) interpreted as the washboard or 
the phase slip frequency . Correspondingly the force / is naturally defined as the work done via sliding 
by one period, that is / =J> /27r/6. Particularly for Electronic Crystals / coincides with the electric field 
strength f = e£ (for Wigner crystals, for ikp CDWs) or with 2e£ (for CDWs, SDWs where one period 
caries the double electronic charge 2e). 



Table 1: Relations between parameters of conventional crystals and CDWs. 



displacements u/b x<^/27r, x = (1, 0, 0) 

velocity, density v = du/dt ^ v ^ —d(p/dt , — Vu dip/dx 

driving and pinning forces / /27r/6 = Qf 

strain, stress \7(p , CV(p 

To quantify and prove the intuitive picture of the Sec. El we consider an isolated local pinning center 
which can be described by a single degree of freedom ipi and monitored by another single one 9i . They 
are the local mismatches of phases 

V'i = </3(xi) - ip , 01 = -Qx; - (p (92) 

relative to the bulk value (p. The latter can be taken to be homogeneous in space, static tp ~ const or 
sliding (p « —vt, within the correlation domain of the collective pinning which by definition contains 
many impurities, see more discussion in the Sec. (HJ. 

Beyond a close vicinity (of a microscopic scale oq) of the pinning center and well within the collective 
pinning domain, we can use the energy functional 



Jv 



2 



(93) 



and typically Vi{(p) = Vi{l — cos{ip)). 

This energy should be minimized over ip{x) at the asymptotic condition ip Cp. By minimizing the 
energy over ^{x) we can get rid of the phase everywhere except at x = Xi. By analogy with electrostatics, 
the "potential" V'(x) originates from the "point charge" V/ = dVi{(p)/d(p, and the elastic energy Wei will 
be the one of a site charged at the potential ip with respect to the infinity: Wei = {K/2C)'4P' where K ^ uq 
(in d = 3) is the "capacitance"^. The case d = 2 is always problematic: now the inverse capacitance 

^Certainly, for the ideal point impurity we face the divergency of elastic deformations which would not allow to determine 
(^(x;): the extremal solution of I93i is divergent at x — > '/^(x) — (p ^ Vj/lxl**""^. As usual, the problem can be regularized 
at a cutoff length ao since microscopically there is a finite width of the pinning site. Actually, there is the microscopic 
coherence length S> no , where the amplitude of the order parameter A ~ ^/C can pass through zero if the phase gradients 
become too high. An even larger length scale appears in quasi Id CDWs: this is the soliton width Is (Figures I16I21I and 
the Section 122) below which the system cannot sustain the shear deformations - see the Section 17.21 
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diverges at the upper limit K^^ ~ \nL where L is a hmiting size. We shall face this effect once again in 
the Section Uni In many respects, the local pinning scheme needs revisions in d = 2 which dimension is 
particular also for the collective pinning, see the previous Sections. 

The elastic regime Wei ~ V'^ is not valid at large deviations where it must give rise to the more 
efficient plastic regime. To see this more clearly, suppose that a very strong pinning center allows for 
the retardation by many periods A^, V' ~ 2TrN, then Wei ~ N'^- The plastic alternative (to emit a pair 
of elementary dislocation loops, the solitons) after each period of retardation - would give the lower 
energy Wpi « 2EsN which grows only as ~ iV rather than as ^ iV^. Actually, for large N a further 
drastic reduction of plastic deformations is possible: as much as from ~ to ~ 7V^/^ In in d = 3 and 
to ^ InA^ in d=2 dimensions, respectively. To see it, remind that the coalescence of dislocation loops 
is allowed provided the total number of embraced chains, that is the total increment/deficiency of the 
crystal periods, is preserved. Then it will be favorable to aggregate the sequence of N emitted elementary 
dislocation loops into the growing single loop embracing N chains which energy is In per unity of its 
length (the perimeter is ~ N^^"^ in d = 3). The expansion of the pair of wide dislocation loops, see the 
Figure [T51 at both sides of the strong impurity will redistribute the retardation by multiple periods along 
the defected line to the retardation by the single period over many lines embraced by the dislocation 
loop. 

We shall return more systematically to the topic of dislocations in the Sec. El 

The above arguments do not tell us yet what is going on within one period of sliding and for inter- 
mediate pinning strengths. Some peculiarities of long range interactions between diverging dislocation 
loops will require to consider carefully also the marginal region between two successive periods. These 
questions will be addressed in the next Section. 



7.1 Basics of metastability. 
7.1.1 Definitions and classification. 

Locally equilibrium states are determined by extremal, over (fix), values of the energy functional at 
presence of only one impurity at the point Xj and for the asymptotic condition tp (p at |x — Xij^cx). 
Here we keep in mind the full scale, generally nonlinear (see examples below) model which is regular at 
small distances and reduces to the elastic model (I93II at a sufficient distance from the impurity. Since the 
ultimately nonlinear pinning energy depends only on the local phase (p{'x.i) + Qx; = il' — 0i (see Eq. (|92|l 
for definitions), it is convenient to keep this value fixed in preliminary stages and optimize for it only in 
the end. Minimization over all other (p{'x.), at given V'(xi) ^nd (p, reduces the energy functional to the 
function. (Below we shall omit the index i addressing only one impurity positioned at Xi = 0.) 

H{^, e) = V{ip -9) + W{ip) (94) 

This variational energy contains the pinning potential V{'ip — 9) (its amplitude we shall call simply V) 
and the energy of deformations W^(V'). It is clear from the above discussion that W{0) = 0, W{2tt) ~ 2Es 
and the next circle starts from this level (the configuration 0' on the figure . so within each period 
W{tp) is the same function with mmW{ip) = W{0) = and ma,xW{tp) = W{27r) = 2Es. The function 
W'{ip), together with —V'{ijj — 9) for several values of 9, is shown on the Figure [T51 Their intersection 
determines extremal values of (|94|l over %p. 

The minima and maxima of the variational energy over -0 at a given 9 determine the branches ^a{9)- 
There are always locally stable states, which can be either absolutely stable, a = 1, or metastable, a = 2, 
the unstable barrier branches arc denoted by a = 3, that is 

n TIT 

anya : ^ = 0; ff(V^(^),e) = Ea{9) (95) 
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Figure 18: Solutions for all branches ipa{0) are obtained by crossing of a thick line (dW/dip) with a thin 
line {—dV{ip — 9) /dip) which are shown for 9 = 37r/4. The doted thick line is drawn for the short range 
model, i.e. without taking into account the long range interaction of dislocation loops. The solid thick 
line shows that these effects lead to the steep fall for the actual W. 
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(96) 



Differentiating Ha over —Cp, that is over 9 along the branch, we obtain the contribution Fa of a given 
impurity to the total pinning force fpin ■ It depends on the instantaneous value of 9 and on the branch 
a = 1, 2 being currently occupied: 



Fa{9) 

27r 



dH 
'd9 



dV 
89 



dW 

dip 



(97) 



Thus Fa is positive/negative for ascending/descending branches, see the Figure [T7I In the Eq. (|97|l 
and henceforth, we normalize the force in an invariant way, as a work performed by displacing over one 
elementary period, 27r in our case of density waves. 

Figure [T7I and Figures [191201 illustrate a typical and more complicated cases. The whole interval of 9 
or some parts of it can be either monostable or bistable, the last case corresponds to the coexistence of 
two locally stable branches: the absolutely stable one with the lower energy Ei and the metastable one 
of a higher energy E2. The same pair of branches can be regrouped also as the ascending branch E-^- for 
which F^{9) > and the descending one E- with F-{9) < 0; they correspond to the retarded and the 
advanced states at the impurity, respectively. Evidently, Ei = min{i?_|_, and E2 = max{i?_|-, _£_}. 
There is a symmetry E-{9) = £'_|_(27r — 9) so that these branches cross at 6* tt, ^^±{11) < tt; also 
Eb{9) = i?_B(27r — 9). In the following we shall assume a certain sign of the overall displacement or its 
velocity, such that branches evolve towards increasing 9 and consider only the important semi-interval 
TT < 9 <2n\ here E+ = E2 and E^ = Ei. 

Notice that the barrier branch has appeared via its absolute energy E3 . But the necessary quantity is 
its increment U with respect to the metastable branch E2 which gives the activation energy for its decay: 



Ui9) ^ E3{9) - E2{9) 



(98) 



This quantity corresponds to the barrier height definition Eb in the collective pinning. 

For strong pinning centers of unrestricted metastability, the two locally stable states coexist over the 
whole period, hence U{9) ^ 0. Typical models show that U{9) is largest at 9 = n and around 2tt (while 
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not exactly at 27r because of special effects of long range interactions, see also the Sec. 19. II below). Then 
there must be the minimum of the activation at some 6 = 9m given by minU{9) = U{9m) = Um (see the 
Figure which plays an important role. 

H 

2E 

s 
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Figure 19: The complete structure of energy branches (thick lines) for pinning centers of highest strengths 
V > ¥2- Uppermost (thin) lines show the barrier energies. Contrarily to the case of the Fig. 1201 all 
branches pass continuously through the whole period. Notice the point 0m where the activation U is 
minimal. The half of the second circle, 6 > 2tt, is also shown. Here, the branches E^,E- are identical 
to those at < < 27r, assuming that the system is totally relaxed. The actual adiabatic continuation 
of the branch is which differs by the presence of two solitons at infinite separation which have 
been created over the previous circle of the branch E^ . The details of the crossover between E+ and E^ 
(dashed circle at the figure) are given on the Figure EH 

For moderate pinning centers of the restricted metastability, the coexistence resides over some intervals 
around tt and 2tt: tt < 9 < 9^ and 9^ < 9 < 27r. At the end points 9,, (or 9^,) the metastable states 
terminate or branch out of barriers, here the second derivative is zero: 

0e, ^e=^(0e): ^=0 = Q ■ U ^ Ve{9, - 9f'^ . (99) 

(see the Appendix Q. The points 6*2, 6*2 appear by splitting off from = tt, 27r at some V > Vi,V > Vi. 
For the short range model (see below) Vi = V\. 

The points 9^, 9e coalesce and then disappear at some higher V = V2 > Vi,Vi (Fig. [T^ which gives 
rise to the point 9„i at V > V2 a,s described above. The full description for V > V2 requires a detailed 
study of dislocations generated near the full period 9 w 27r. We shall postpone the analysis of all special 
effects related to diverging dislocation loops or solitons till the Sec. 1^1 

All together our qualitative arguments lead to the structure of energy branches shown on the Figures 
I19I20I There may be a more complicated picture of termination points for less local pinning centers, e.g. 
the Figure and even a more complex hierarchy for the collective pinning regime. At present, this 
feature of the energy landscape stays beyond the scaling theory of the collective pinning. 

7.2 Models. 

The above analysis exploited only the most general properties of the variational energy H: the periodicity 
of the pinning component Vltjj — 9), the monotonous energy of deformations W^ip) with the minimum 
at tp = 0, the maximum at ?/' = 27r and with the inflection point ipm in-between; and complemented by 
universal long-range effects of distant dislocation loops (Section^. Now we shall illustrate these features 
on model examples. 
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Figure 20: The complete structure of energy branches Ea{0) for pinning centers of different strengths: 
weak {V < Vi, lowest thick curve) and intermediate {Vi < V < V2, other thick curves). Uppermost 
(thin) lines show the barrier energies. For intermediate strengths, notice the disconnected region of 
higher energy branches around 9 = 2-k, in addition to continuously accessible ones corresponding to the 
Figure El The termination points of the two regions 9e and 9e coalesce at V = V2 giving rise to the 
structure of the regime V > V2 shown on the Figure [T^ 



7.2.1 Elastic model 

This model takes into account the periodicity of the pinning potential but neglects the topological char- 
acter of plastic deformations derived from the same lattice periodicity. As a model of metastable states 
it was already considered in |18| . For the bare pinning energy it is always natural to choose the point 
impurity pinning potential V{ip) ~ V{1 — costp). The deformational part W{Tp) is the energy for the 
distribution ip{x) optimized at the condition that 1^9(0) ~ ip + and ^p{oo) = (p. To calculate this energy 
we consider a small sphere 5) of the radius ao around the position of the impurity on which we assume the 
phase is (/5 = ^ + We have then to solve the Poisson equation V^<y9 = with the boundary conditions 
on the collective pinning domain V far away from the impurity {Ln ^ oq) and on S. The solution and 
the energy are, in d = 3, 

0(x) « <^ + , """^ , , WitP) = 2TraoC^^ , k = V/iinaoC) (100) 
|x - x^l 

Depending on the relative impurity strength there is only one {k < 1) or more {k > 1) solutions of Ea. (|95|l 
for the energy branches 

dH/dtp = , Ksin(V' - 6*) + V 0. 

The condition k = 1 identifies Vi = AttgoC and the results of the above general treatment follow corre- 
spondingly, except for the region V > ¥2- The latter requires that the periodicity to be fully taken into 
account: W ^ tp"^ does not show the inflection point tpm- 
Detailed calculations for this model can be found in |108j . 
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7.2.2 Solitons in quasi Id system: the short range model. 

Wc consider a quasi ID system of interacting CDW-chains with an impurity at the chain n = at the 
position a; = 0; the Hamiltonian is 




(101) 



Here the first term in square brackets is the on-chain elasticity and the second term is the interchain 
couphng which is reduced to the shear elasticity ^ Cj_(V_l(^)^ when perturbations are small, C_l = 
Tlim^mn- The 2tt pcriodicity of the pinning energy allows to skip the 27r quanta in tp^ to optimize the 
total energy which already can originate local metastable states as we have discussed above. Moreover, 
the 27r periodicity of the regular energy in (|101|l allows for topological defects, the solitons. For the soliton 
centered at the position X at the chain rt = 0, the phase profile ipg{x~ X) describes stretching/dilatation 
by one period along the chain ?i = relative to surrounding chains, see the Figure |^ It is distributed 
over the length Is ^ \/(^\\ I C'± and costs the energy Eg ^ ^C||Cj_, the last defines their equilibrium 
concentration ~ exp{—Es/T). 




*^ soliton 
adjustment 



Figure 21: Extinction of a point impurity pinning at presence of the 27r-soliton. The phase profile 
ip{x — Xs) can be adapted (the vertical arrow) to the phase mismatch at the impurity position Xs by the 
adjustment (the horizontal arrow) of the soliton position Xs- 

While single solitons can be created only by phase slips, their pairs are non topological configurations 
which can be continuously developed by driving 9. Hence we are looking for extremal values of (|101|l with 
trivial boundary conditions ip{±oo) = 0. They can be visualized (see the Figure ITC)) as a combination of 
two pieces of ±2tt solitons at positions ztX: (ps{—X) — (fisi^) ~ from which one concludes the relation 
X = X{ip). We can specify the shapes (ps{x) within a short range model iSSj suggesting that in (|101|1 
only the central chain n = (passing through the impurity) is perturbed while its Z 3> 1 neighbors stay 
at ip„-io = <y5 homogeneously. Then the energy functional is simplified as 



dx 



2^ll(^') 



f\2 



C± C0S(ly9) 



V cos{ip ~ e)6{x) 



(102) 



Its extremum is the function ipa{x — X) — ips{x + X) where fs{x) is the standard saddle point solution 
for the Sine-Gordon soliton. The sequence of plots of (p{x) for different 6 is shown in the Figure ITCl The 
energy is W{'tlj) = Es{l — cos(?/;/2)); in accordance with our general expectations, over one period W{tp) 
changes monotonously within = W{0) < W{'iIj) < W{2n) = 2Es, W{tp) ~ tp^ at ip ^ 0. Now we can 
identify the threshold values as Vi = Eg/ A, V2 = Es/2, ipm ~ ■n 12. Typical solutions of the equation (|95|l 
for this model are shown on the Figure |221 

The short range model (|102|l already contains most important features necessary in applications of the 
local pinning. It fails only for high velocity regimes of very strong impurities when the two well formed 
solitons diverge at ±X — s- ±00. In this regime their interaction with each other and with the impurity 
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Figure 22: Solutions tpa{d) for the short range model shown at < < V2. They correspond to the 
plots of the Fig. 1201 The double period for — 27r < 1/; < 27r allows to see both retarded (medium lines) 
and advanced (thick lines) configurations; thin lines correspond to barriers. 

penetrates very efficiently, as a power law, via the elastic deformations of the whole media. For the short 
range model, where the surrounding chains n ^ are frozen, the perturbations fall off exponentially 
as exp{—X/ls) which gives incorrectly the analytic phase dependence of the energy W{'ip) — 2Es sa 
£'s(V' ^ 27r)^ at ■(/; ^ 2ti. Contrarily, the true power law for long-range elastic interactions results (for 
d = 3) in W{'ip) — 2Eg ^ Es{ijj — 27r)^/^ which leads to particular instabilities. This effect relies upon the 
view of solitons as nucleus dislocation loops, and we postpone its analysis till the Sec. 19. II which collects 
all information on special contributions of dislocations. 

8 Kinetics and relaxation, v-f characteristics. 

Results of the previous chapter provide the basis for the picture of the local/strong pinning. In this 
chapter, we shall apply this picture of local metastable states to kinetics, and finally we shall describe 
nonlinear v — f characteristics and to the linear response function x(w). 

Before going on, it is appropriate here to rectify our definition of the local or strong pinning. When 
these notions were introduced |18II2()| . the strong pinning case implied primarily that the local adaptation 
of the elastic media follows closely the minimum of the impurity pinning potential. Our definition compiles 
with this tradition but generalizes it to the multivalued case when the minimum is allowed also to be the 
metastable one. But the same time we disagree with a rather common view that for the strong pinning 
the correlation volume is of the order of the one per impurity. At least at d > 2, we see that however 
strong the impurity potential is, the deformation falls off with distance as a power law. At d > 4 the 
perturbation would be levelled out completely while in d = 3 it will contribute to long range fluctuation 
of the collective pinning for which there is no difference with respect to the impurity strength. The 
next ambiguity is to identify the strong pinning via the linear dependence of the critical field fc on the 
impurity concentration n^. This definition is also traced back to the old epoch when the importance 
of time scales was not appreciated yet. Today it is clear that, because of limited heights of potential 
barriers, the local/strong pinning describes relatively high velocities or frequencies, or respectively low 
temperatures, where its contribution is indeed linear in rii. 
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8.1 Kinetic equation. 

Wc consider now the kinetics in an ensemble of impurities possessing local mctastable states. With the 
exception of particular regimes of strongly divergent dislocations, these states are formed locally, at a 
distance shorter than the mean impurity spacing Li which, in its turn, is much smaller than the pinning 
length Lp. Hence we can define a reference phase for a volume D staying well within the collective pinning 
correlation volume but still extending over the large number of impurities: 



(^(x, t) = V-^ / rfx(^(x, t) ; Lf < P < Lf' 
Jv 



IV 

Here we can neglect the dependence on x, so that (p « (p{t). For the same reason wc can neglect the 
direct contribution ~ /oq of the driving force to the energies Ea- Moreover, we can separate the local 
pinning force /pm from impurities and the regular viscosity fvisc ^ v/j coming from the phenomenological 
damping (|39|l : / = fvisc + fpin Here the time delay related to viscosity St ^ clq/iC is small in compare 
to . For the collective pinning this problem is more complicated because the microscopic scale ao is 
enlarged to an intermediate one Lp (see Eq. (I65f) and the related text afterwards). In the following we 
shall speak only about the pinning part fpin of / implying that fvisc can be added in the end. Actually 
even that is not necessary: we shall see below that the linear damping ~ u is generated by the pinning 
itself which result is confirmed by experiments, see the Sec. [TSl 

As in statics, each impurity can be characterized by the positionally random phase 6i = — Qx; — <^(t) 
which now evolves in time following the moving reference phase (pit). The single monitoring parame- 
ter describes both the time evolution and the distribution over impurity positions; the average over 
randomness, {■■■)r becomes 

<■■■■>- 

For each bistable impurity (see the Sec. [7|and the Figures [17119120(1 . the state occupies instantaneously 
one of the two branches "±" with energies E±{9): ascending F+ > and descending F_ < where F± = 
2TTdE±/d9 are the local pinning forces. The upper and lower energy branches arc E2 = max{i?_, 
and El = mm{E-, E+}. There is also the barrier branch Es^O), E3 > E2 > Ei. The branches E^, E2 
can terminate at points ^e, Stt — ^e, where both the upper metastable branch E2 and the barrier £^3 
split out: i?3 = E2 at 9e, hence at 6'e there is no activation energy, U ~ E^ — E2 ai 9 ^ 9^: for 
the decay of the metastable state E2 to the stable one Ei. The barrier activation disappears at 6'e as 
U ~ Ve(9e — 9Y with = 3/2. The branches cross at 6* = tt and we neglect the repulsion between Ex and 
E2 due to quantum tunnelling, see the Section [TUI In the course of the density wave motion 9 = 9{t), the 
distribution of occupation numbers n± = {n_|_,n_} for branches ± 

n± = ^ilTn); n^ni9,t) (103) 

obeys the kinetic equations (see more in the Appendix Id|) 

dn rieo — n d d d d9 , /S.E ^ ^ 

v; n,q = tanh — , l\E ^E+ - E_ (104) 



dt T ' dt dt 89 dt ' ^ 2T 

where Ueq is the value of n in thermal equilibrium. 

Here and mostly below, we imply an internal relaxation which is due to passing over the local barriers. 
Its rate t{9) does not depend on the velocity 9 but is an essential function of the position 9: t ^ 
To exp(/7/T) where U = U{9) and Tq^ is an attempt rate. Later, in the Sec. El we shall discuss also the 
external relaxation which is due to the mean free path of diverging pairs of dislocation loops or solitons. 
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Finally, the pinning force averaged over both the sliding period and the initial conditions is given as 



where Oieft, bright = Stt — are the bistability limits. In the last form of /, the first term niAE\g^^^^ = 
fmax gives the energy dissipation by the ultimate falling from the termination point, while the second 
term erases this value fmax due to occasions of earlier fallings down which arc more frequent at lower v. 



8.2 Stationary motion. 

Consider the stationary process when the density wave moves with a constant phase velocity v = ^(f = 
9 = const, then dn/dt = 0. The solution of the kinetic equation H1U4|) is trivial, see Eqs. 1)1401 114l|l 
in the Appendix IdI being even simpler at low T ^ AE{6): For AE{9) 3> T, the pinning force can be 
written as a weighted distribution of instantaneous forces: 

dOFie)o.v I - / £^ 1 , = = (106) 




or of energies AE dissipated via falling from the metastablc to stable branches: 

d9i \ 



fmax CXp 




vt{9i ) J 
d9 ^ / d9i 



2.n.l -^^AEi9)e.J-l 



vt{9) ' ' ' \ J. vt{9i) J ' 
fmax = 2nn,AE{9 

max 

Here fmax is the largest value of the pinning force to which it saturates at high w; 9max = ^e, 27r (depending 
on the impurity strength) is the most distant point reachable by the metastable branch. 

The Eqs. (|106|l and H1U7|) are suitable for calculations at small and large velocities correspondingly. 
After simple calculations presented in the Appendix^ we arrive at the results shown schematically at 
ihev- f plot of the Fig. ESI 

1. Small velocities v ^ t~^: 
where = '''('''') ^ cxp^Utt/T), U-„ — U(tt) is the maximal relaxation time in the region of the branch 
crossing point 9 ^ tt. The main contribution comes from the close vicinity of tt: 9 k tt + 59 where 
59 ~ vt{9). We can distinguish between two subregimes. 

la. Very small velocities v <^ z;^ = {T / Ftj)t^'^ <C t^^, — F(tt) : 
The decay happens as soon as the branch becomes metastable in a vicinity of tt, even before the 9 
dependence of t is seen. The life time interval is 59 ^ wr^, hence the law 1)106(1 yields 



dAE 

f = nriiVT^ 



d9 



(108) 



which emulates the phenomenological viscosity. This is the regime of the linear collective conductivity 
c = ^// = const. It shows an activated behavior via which can emulate the normal conductivity. 
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Figure 23: Schematic plot of the f{v) dependence showing several regimes of relaxation, see the SecIHl 
for explanations, here for shortness fmax — > fm- The zoomed vicinity of small w = and f m fc should 
recover the collective pinning sliding regime with an opposite curvature of f{v) and finally, at very small 
V of the collective creep regime, the curvature will change once again (according to the Figure [71). Note 
that the viscous force is not included in this figure; it would simply give an inclination to dashed line of 
the asymptotic regime. 



lb. Moderately small velocities v-^ <^v<^t^^: 
Deviations 60 are still localized around tt but they may be already large enough to see the decrease 
of the barrier height: U ^ U„ - F„S9/2tt, t = cxp[-(5e'i^^/(27rT)] , F^/T > 1 Wc have exploited 
the fact that the branch Eb{0) has a minimum at tt so that the linear dispersion of U{6) is given by 
dE2{0)/d9 ~ Ftt. The condition vt{9) ^ 66 is fulfilled at 69Ft^/T ~ ln[i;r7r] and finally the dependence 
f{v) or v{f) become 

/ n.iT\n[vT^] ; v ~ t~'^ expif/n^T) 
Convenient interpolation formulas for cases la,b are 

/«27rTnan(l + .|^) ; . = |^ (exp ^ - 1 

but the integral representation p47|l is more precise. 

The physics of / ~ w regime is given by the high probability to stay on the metastable branch during 
a small displacement 69 Tt^v. The physics of / Inw regime is that at higher i; a wider region of 69 is 
explored and the metastable branch starts to feci the decrease of the barrier (long ahead there is either 
the termination point 9f, or the minimal barrier point O^m even if unreachable yet at these moderate i;). 

2a. High velocities v 3> t^t^^.' restricted metastability Vi < V < V2: 
At higher velocities, wr^ ^ 1, the points distant from 9 ^ n can be explored, and especially 9 9^ 
becomes important. The motion along the branch starts to reach a close vicinity of 9^ at high enough 
Ve ^ V ^ where, see ()149|l . 

TX'/" 1 3 



, , ^ ^-^V,^ cnst (109) 

Ve/ To Z 

and the coefficient T4 was defined in (|132|l . Within the limits 3> w 2> Wtt: the diminishing barrier 
U « Ve{69Y^^ still stays at U ^ T so that the decay of the metastable branch E+ is still activated. 
Finally we obtain from H107|) 



f = fmax - CeUiFe — In — ] V Ve CXp 

Ve V 



Ve f fmax f 



T \ CeJl.Fe 



Ce ~ 1 (110) 
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Realistically, this regime can be found only at very low T: the crossover velocity must drop well below 
the microscopically high values of the attempt rate Tq"^. 

2b. High velocities wr^ S> 1.' unrestricted metastability V > V2. 

For very strong impurities the metastability is maintained over the whole period. Both the metastable 
branches and the barrier branch stretch over all 9 and t 7^ everywhere. The important role is played 
now by the point 9m where the barrier activation is minimal: min?7 = U{9,n) = Um- It gives rise to the 
minimal (over the whole branch) relaxation time t„i ^ cxp(Um/T), see the Figure [TUl In the vicinity of 
9m we can write U{9) = + {9 — 9m)'^B/2), B ^ 1. There are two different regimes described below. 

3a. High velocity range v ^ Vm = Tni' ~ max{r~^}: 

Now the 1 /v expansion works in H107|) and we find from p5()(l 



n / T 

f = frnax ' COnst^ J —— [/^E{2ti) - AE{9m)). fmax = 27rn,;A^(27r) (111) 
VTm V BUm 

The asymptotic force fmax is the energy to create the pair of dislocation loops (solitons), it is determined 
only by the final point 9 = 2tt. It is approached by the law / — fmax ~ — l/w which reminds, at first 
sight, the collective pinning corrections for high velocity, but the sign is opposite! 

3b. Moderately high velocities u,„ 3> w S> t^^: 

Now V is high enough to reach the point 0m, but not that high yet as to bypass it easily - still VTm ^ 1, 
hence the point 9m will provide the major relaxation. This effects will be particularly pronounced near 
the threshold V2 when Um, and consequently Tm, become much smaller with respect to their typical 
values over the branch. We can easily obtain (up to numerical coefficients) 



/ = 2nn, {AE{9m) + (AS(2^) - A^(0,„))P} , P ~ exp ( — ] 

y V BUm VTm J 

This formula tells us that the main force is provided by the part of the branch between and 9m or, 
equivalently, by the energy released from the relaxation at 9m - similar to termination points for the case 
lb. The second contribution ~ P comes from the remaining part of the branch, between 9m and 27r, 
but the penetration probability P to this part through the "hot point" 9m is exponentially small in 1/v. 
Nevertheless this small probability is responsible for the irreversibility and memory effects. Indeed P is 
a probability to create dislocation loops (solitons) which are long living plastic deformations. 

Summary: The most important cases of these regimes are shown on the Figure !^ and will be discussed 
again in comparison to experiments in the Sec. ^1 We have skipped from consideration the most limiting 
cases: highest v for strongest U will be considered in the Sec. devoted to effects of dislocations; lowest 
V (the law on the Figure ESI) will be considered in the Sec. ITTl devoted to the ensemble averaging. 



8.3 Linear response. 

The standard response function x = ^'P/^f is measured in CDWs as the dielectric susceptibility e ^ x- 
Within the collective pinning picture it can be found with the help of the same kinetic equation as 
described in the Appendix ^ Here we shall follow a more intuitive and transparent approach first 
considered in Ref. [21]. Consider the reaction of local bistable states to a weak perturbation for low 
w or at large t. In equilibrium, the impurities occupy the lowest branch Ei (for a moment we neglect 
the effects of the thermal population). For small variations 69 the main contribution comes from the 
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degeneracy point 9 ~ t: where the two stable branches cross each other. Impose a small homogeneous 
shift of the relative impurity positions 69 at t > 0. In the whole ensemble of impurities, those ones 
become metastable which have been occupying the interval of stable positions tt — 59 < 9 < tt while the 
stable positions at n < 9 < n + 59 become empty, as shown on the Figure El The subsequent evolution 
follows the relaxation towards the thermal equilibrium. The imbalance of forces gives rise to the inverse 
response x~^. Going from the real time to the Fourier representation we find 

5f{t) = 59n,F^exp{-t/T^), 9^^-, Xj^ = ^ = ^ ^"'f ■ (112) 

At high oj ^ ^ Xu> saturates at its maximal, real value maxx ^ (riiFTr) ^ ■ The small u; limit of x(^) 
corresponds completely to the small velocity limit of the v — f law (|146|l . At ojt^ ^ 1, the system reaches 
the thermal equilibrium and xZ^ disappears ~ iu giving only a contribution riiFT^TT^ to the damping 
parameter 7"^ in full agreement with the f — v results. 



9 Generation of dislocations at high velocities. 

By now we exploited mostly the general properties of metastable branches: the existence of points of 
levels degeneracy AE ~ 0, of the barrier termination U = or of its minimum. Topological defects 
were implied to exist in the background providing peculiar reasons for the metastability. Quantitative 
results were derived for a general position of 9 when dislocation loops have not emerged yet as distinct 
(and distant) entities. A more elaborated analysis is required near the final point 2tt of unrestrictedly 
metastable branches accessible at high v. Here the metastable configurations are formed explicitly by 
divergent pairs of solitons, more generally dislocation loops, centered around ±X - see Figures 1 151 and 
[T51 Their interaction will modify both the structure of energy branches and the related kinetics. 

First of all, we add here a few technical notes necessary to work with dislocations. For details see 
HSllHSlEni in general and [^2 EI specially for CDWs ^. In CDWs particularly, dislocations have all 
properties of conventional vortices in planar magnets or superfluids (with an exception for a special 
conservation law for the total area embraced by the loops which distinguishes the climb from the glide) . 
Even more pragmatically, we can invoke a common wisdom of the magnetostatics considering dislocations 
as currents, the strain as the magnetic field Wip ^ H and the stress as the induction CV<p 4^ B. The 
signs are however different: antiparallel dislocation lines (dislocation lines with opposite polarities) attract 
each other. For our typical case of dislocation loops lying within the plane (y, z), i.e. perpendicular to the 
sliding axis x, we arrive at the following prescriptions: In the inhomogeneous field (/j(x) created by other 
sources, e.g. another dislocation loop or the impurity, the glide force in x direction applied to the unit 
length of the dislocation loop is ~ CW±ip. The energy per unit area of the dislocation loop is ~ Cdx^p- 
The dislocation loop self-energy has the standard vortex form: at large R, Edl{R) ~ Ci?ln(i?/ao) in 
d = 3 or ln(i?/ao) in d = 2, respectively. At the minimal i? ~ oq the dislocation loop is interpreted as 
the nucleus dislocation loop embracing just one chain - the soliton, E£>l{R) Eg. The phase distortion 
by the dislocation loop itself at a given point is the half of the steric angle (the angle in d = 2) 

at which one views the dislocation loop form this point. This angle evolves from to 27r along the path 
crossing the dislocation loop, these asymptotic values are approached as ±5ild ~ {R/ X)"^^^ . 



9.1 Effects of dislocations upon metastable states. 

Consider firstly the long range instabilities near ?/;, ^ = Stt corresponding to the divergent pair of solitons. 
From large distances, the dislocation loops interact with each other and with the impurity via long 

®For more complicated techniques of working with ensembles of dislocations see 11091 . 
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range elastic forces. By definition, the presence of dislocation loops at points ±X displaces the phase 
in between, at x = 0, by 2t: with a deficiency —Sip = 2tt — ipiO). The latter is given by the sum of 
(steric) angles ftd ~ {R/XY~^ (at large X/R) of their view from the point x = 0, see the Figure [T51 
The energy W is equal to the energy of two dislocation loops 2E£)l{R), also taking into account their 
mutual interaction at finite X. The attractive potential 6W of the loops is given by the longitudinal stress 
Cdx^p ~ —{C / R)[R/2X)'^ produced by one loop over the area ~ R<^-^ of another one at the distance 
2X: 5W CR'^-'^{R/2XY. Finally, we exclude X in favor of -0 to arrive at 

5iP - --{R/Xf-^ , SWii^iX)) - ^CR'^~^iR/2Xf - -CR'^-\~6i:f^^'^~^^ (113) 

We notice that, while the force disappears as W' ~ (—5ipy^^'^~^^ 0, its derivative diverges W" — *■ oo 
giving rise to the branch instability. Namely, even for arbitrary large V ^ W we shall meet the condition 
W" = —V", hence there is always a solution of the last equation in l|99|l at some ijj* < 2ttJ A more 
elaborated analysis given in the Appendix |0 shows that this value of ?/; is reached at 6' = 6'* > 2tt, that is 
already in the next circle of the mean sliding. According to the first relation in the Eq. H113() . there is a 
distance between dislocation loops Xg associated to the phase deficiency V'e — 2tt. We see that instead of 
diverging {X ^ oo, which requires for ■(/i 27r) at the end of the period 9 211. the pair of dislocation 
loops looses its stability at some critical distance X* corresponding to V'e < Stt of the branch E'^. Then 
it falls to the higher distance X* => X^ corresponding to tpe > 2tt of the branch and continues 
the new circle along this branch. As a result, the region encircled by the dashed line on the Figure [THI 
acquires the structure shown at the magnified picture of the Figure !^ We shall call 6*^ the overshooting 
termination point of the overshooting part E*^ of the branch E+ penetrating into the next period 6 > 2tt. 
The appearance of new termination points brinks some similarity to the vicinity of the tt point at y > l^i . 
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Figure 24: Special effects near the crossover between two successive periods which are originated by 
interactions of distant dislocation loops. E*^ is the overshooting part of the branch E+. The branch E+ 
differs from the lowest branch E+ by presence of two solitons at ±oo. 

The possibility for relaxation of E'^ to E^ is not unique. In principle, there is always an option to 
fall down directly to the lowest branch E+ (the same 9 but if) being close to rather than to 27r), thus 
releasing the energy k, 2Es. But it requires for the annihilation of the distant pair of dislocation loops 
which involves a large barrier both in energy and configuration. There is one more option: to fall to a 

^In CDW materials there are also long range Coulomb forces which are not screened at low T. They affect drastically 
the energetics of dislocations 1941 which leads to an even more singular law W^'ip) — 2Ea ~ Eg\ip ~ 2tt\. 
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close, in energy, descending branch (see the Figure E^ . But this is the strongly advanced branch 
(all descending branches are advanced configurations) corresponding to the pair of dislocation loops with 
opposite polarities, so that this transition would require to switch phases from w 27r to w —2-k along the 
whole interval {—X, X) which is hardly possible. 



9.2 Kinetic effects of diverging dislocation loops. 

The complete kinetics of these states is complicated for several reasons, one of them is a larger number 
of branches involved. A simplification comes from the high velocity condition to reach this regime: 
V 3> \/t{6) at all 6. Here r is the relaxation time for dropping from the term Ej^ to which we always 
kept in mind before. Neglecting this basic relaxation in a small vicinity of 27r, we can concentrate on the 
short relaxation time r* ^ r to fall from E^ to Ej^. Now on top of the law / — fmax ^ —riiF/vT (the 
case 3a of the Sec. 18. 2f) we shall see the sequence of regimes similar la, lb, 2a (also from the Sec. I8.2|l 
for the case of the restricted bistability, but with much smaller r* instead of r. This new contribution 
to f{v) falls off slower (hence finally winning) than ^ \/v but its emergence is delayed because the force 
is reduced to the smaller value F*, as given by the inclination of the overshooting branch E'^ (see the 
Figure E^ . We shall notice traces of this regime in applications to CDWs in Sec. 1131 

Actually, the relaxation time approach may not be applicable any more. For well separated solitons, 
the extrinsic mechanisms of relaxation enter the game: annihilation of solitons and antisolitons produced 
by neighboring pinning centers along the chain, aggregation of solitons into growing dislocation loops, 
disappearance via phase slips. Now the final rate is determined by the soliton distance X in comparison 
to its collision mean free path A rather than by the time in compare to r as it was for intrinsic processes. 
Phenomenologically, A is included to the starting kinetic equation by the following substitution which is 
noticeably different from 1)11)4(1 : 

d d 1 dX{%lj(e)) 
dd~' de^ A d0 ' 

where X{ip) is the distance associated to the retardation ip{9) taken along the term E+. It results in the 
following modification of Eqs. HlUt)llU7|) . as well as of (|141|) : 

de fdOdX X(9) , , 

- — I - ' ' (114) 



^ vT{ei) J X de A 

Remind now(see Eq. H113|) and the text above it) that the diverging dislocation loops of the radius 
R leave in between the phase retardation approaching 2tt as (at d=3) 6ip = il'iX) — 2tt ^ — i?^/Ar^; and 
the same long range elastic deformations provide their attraction with the potential ~ —CR^/X^. Well 
before the overshooting instability develops, that is at 59 < 0, we have (see the Appendix lO ^ 5) 5iIj ^ 59 

as it should be for a very strong impurity. Then 5n ^n-2Es Es{-59)^/'^, hence F - Es{~59)'^/'^ , 

and we arrive at 

i? R 

X ^ — . : s ' 



-59 ' X^/^ 

The condition s 1 in (|114|l defines the characteristic 59 ^ R/X and finally we obtain the force correction 

5f,na. = -n^ J 2F^ ^ -n,Es{-5ef'^ ~ -n,E, ' (115) 

The last formula relates the high v asymptotics of the pinning force and the mean free path A of dislocation 
loops. 



50 



10 Quantum effects. 



The strong pinning picture gives also an access to quantum effects which become important at very low 
temperatures, when the thermal activation is not efficient. The quantum creep became the subject of 
experimental studies since jllOj which work has attracted a substantial attention in theory, see e.g. 
The existing schemes concentrate upon the quantum nucleation of CDW advances in regions free from 
impurities. But this approach cannot tell us how the pinning is released and the motion as a whole is 
initiated. For this goal, it is necessary to consider the quantum decay of mctastablc states at the pinning 
centers. 

Phenomcnologically. the dynamics is introduced via the kinetic energy {I /2)ip'^ with / being the 
"moment of inertia" associated with the "angle" V'- Typically, / ~ ljq'^ where is a microscopic scale 
for frequency of local quantum vibrations, usually it is associated with the attempt frequency wq '''o^- 
Then the quantum Hamiltonian is 

H = H{tjj, 6) = Pl/2I + H{ij, 9); P^ = -ihd/dip 

where is the momentum conjugated to "0 and 9{t) is considered as the time dependent parameter. 
The quantum interference is efficient only near the branch crossing point 9 = tt. Then the degeneracy 
will be lifted by splitting the levels Ei and E2 which opens the quantum gap Sq between them as shown 
on the Figure 1^ Within the normal dynamics of Eq. (|10|1 . the gap is Sq ^ exp{— cnst\/lU /!%). But the 
emittance of phase phonons to the bath will have a pronounced dissipative effect upon the tunnelling as 
we shall describe in the end. 

Hi 




Figure 25: Pinning extinction by the quantum tunnelling between branches E±. The gap Sq opens 
between classical branches Ei and E2 which were degenerate at 9 = tt (compare to the Figure [T7|) . 
Adiabatically, the system follows (the arc arrow) the lowest branch Ef, thus giving the zero force in 
average over one period. Only nonadiabatic transitions (the diagonal arrow) from Ef to allow to 
reach the metastable branch to gain a net pinning force. 

Working within the normal dynamics, we arrive at the standard Landau-Zener problem of tunnelling 
due to slow time dependence of the Hamiltonian (via 9{t) in our case). The standard notion is that if v, to 
are negligibly small, then the system will follow adiabatically the exact quantum branch of the lowest 
energy Ef which is a mixture of classical branches Ei^2 

El ^E2 + Ei- ^{E2-Eif + Sl (116) 

and the average force is zero within an exponential accuracy. Namely, at low v <^ Sq the force will be 
determined by the small probability of the nonadiabatic transition to the upper metastable branch E2'. 

f ~ exp ( —] ■ —^Sq 

\ VTqJ Tq 
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Oppositely, the tunnelling between bare branches £'1.2 is suppressed for large velocities or frequencies 
of 9{t) (v,uj 3> Sq), and the average force is given by the classical picture we have described above in 
terms of the thermal relaxation time. In kinetics, the effects of quantum and thermal fluctuations seem 
to be similar and in our simplified picture the inverse times are additive t^^ = + t^^. 

The phenomenologically introduced inertial dynamics of the Eq. HlOfl may not be valid. Actually, the 
dynamics of the variable ip becomes retarded and dissipative because of the nonlocal contribution from 
the whole field (y9(x, t). Its action (on the imaginary time axis) is given by 

Sbuik['p{^,t)] I J d^^dt^u-^ {dt^f + (V^)2} 

where u is the phason velocity and C is the static elastic modulus. Now we can repeat the prescription 
of the Sec. 17.11 to integrate out iy9 at all x 7^ keeping fixed the value 1^9(0, t) = (p(i) = -0 + 6*. In the 
Fourier representation we have 

Sk^nW{t)] = ^ ^ |(^^|' , = I^U^ (117) 

u 

where /^^ is the frequency dependent generalization of the constant /. while the kernel is given by 
C f dk 1 f dk ( 1 / 1 1 



J (27r)''k2+cj2/u2 J (27r)'' \k^ k^+uj^/u^ 



(118) 



Mainly we shall address the case d = 3. Here we should use the second form in the Eq. H118|l where 
the first term in brackets {} gives, after the regularization at high k ~ ^, just the elastic contribution 
H100|l we considered in the Sec. 17.21 its scale is Kq ~ Cqq. The second term, regular at high k and 
hence model independent, gives the contribution ~ \oj\/u and finally we obtain ^ CaqIwI/u, while the 
regular contribution ~ appears only as the next order in perturbation. The form IjllTf) with ~ 
is typical for the dissipative quantum mechanics |112| . which route we shall follow below. In the time 
representation, the gives rise to the kinetic action 

s,Mtit] ^ nr fdt.dt, ( m^^\" , r = ^ » 1 (119) 







ti — 12 J hu 



It is logarithmically divergent in the tunnelling time t if the tunnelling trajectory acquires a final increment 
(between tt and — tt in our case). The total action can be written schematically as 

/ u 

S{t) = hr\n{tLJu) H \-Ut , uou^ — 

t ao 

Here we have included also the regular inertia moment / for which there are always some local sources. 
The tunnelling level splitting is Sg ~ exp {—Smin/f>) where S„iin = miniS'(i). 

At intermediate 1 <C F <C V lU / h, we arrive at the usual WKB result Smin ^ \^IU, but with an 
essential preexponential suppression of tunneling: 

' u V^'^ ( VTu\ 



cxp —cnst 

At higher F S> VTU /ft, the tunneling suppression is more drastic: 



Smin ~ F In ( I , 5q 



r 
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The last condition imposes the constraint upon the value of the pinning potential which must be com- 
patible with the metastability condition. For typical models we find ujo ^ to be necessary. 

In d = 2, the first form of K~-^ in (|118|l should be used. Now the whole integral is diverging at small k 
yielding a universal result. We obtain an even slower frequency dependence ~ C/ In {u/\aoijj\) that is 
K{t) ^ l/{t\-n^ t) instead of as in H119(l . The logarithmic divergence is the same we have noticed for 
the static problem. We see once again, remind the Section FT^ that in d = 2 short and long range effects 
cannot be separated, whatever they concern: the interference of the collective and the local pinning, or 
the local dynamics and the one related to emittance of sound to the bulk. 

Apparently further studies are necessary. Already now we can understand the quantum creep as 
the tunnelling between retarded and advanced configurations at the moment when they become almost 
degenerate. The process is strongly affected by emitting sound excitations which drive it to be dissipative 
even at T = 0. 



11 Ensemble averaging of pinning forces. 

Above, in studies of both f — v and Xj we have assumed the simple exponential relaxation at identical 
pinning centers. In reality, there may be either a broad distribution of impurities strengths or a tail in 
addition to the peak at the value for a typical pinning center. Effects of distributions can be important in 
applications and they are particularly necessary to build a bridge to the collective pinning regime where 
the broad distribution is the basic ingredient. We shall concentrate on the most pronounced effects at 
lowest V and w compatible with the local pinning picture. 

For a distribution of barriers VuiU), the distribution of r is T',- = VudU/dr = VjjT/t and we shall 
consider two limiting cases. Firstly, the model with the exponential distribution of barriers corresponds 
to microscopic fluctuations of e.g. a distance between the CDW chains and impurities: 

Vu - C/o"' exp(-C//C/o) , Vr ~ (r/C/oro)(To/T)i+^/^°. (120) 

Similar effects appear for Poisson and Gaussian distributions. Secondly, we can try also the scaling 
distributions which appear intrinsically within the collective pinning regime, now Vt ^ (rln" r)~^ where 
the index a > 1 (this condition is necessary for convergence of the normalizing integral) depends on the 
dimension d and the critical index x)- Naturally, the distributions are normalizable, but we also notice 
that in all cases their first moment, which is the mean value of t, is divergent: 



J drVrir) = 1 but (r) = J dTVr{T)T 



OO. 



Remind now that the low ui, v asymptotic for both xZ^ sjad f{v) are linear in r (xZ^ ~ ^^'^ fi'^) ^ '^''') 
and then saturate or change to a slow growth at higher crossover values wr, vt I. Hence, their averages 
will be divergent within the ~ r regime and saturate at the crossovers. E.g. for the response function 
H112|l we obtain 



/iujT 



QxZ^ - n^F^Lo^Vrito ^) ^n,F^Vu{THu ^)) (121) 

Thus, the low uj tail of the imaginary part 3%"^ gives the direct access to the distribution of potential 
barriers. A similar result is obtained for the real part 3?xJ^ which is given by the integral of r^, the second 
moment. The interpretation is that at relatively low uj (still within the local pinning domain) or high 
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T, only those long living states contribute which are due to rear occasions of large impurity potentials, 
hence large barriers U ^ Tlnoj^^. Then ~ VuiThi l/uj). E.g. for the exponential distribution (|12()() 
we find xZ^ ~ 

The same procedure can be applied to the f — v dependence. The law 1)1 08(1 was derived for a typical 
impurity. At the lowest v, we may still find some rare regions - may be clusters of host imperfections or 
particularly strong impurities - where barriers U are high enough so that vt(U) ~ 1 still holds, hence 
we are looking for U ^ Tin (vq/v) and the pinning force is given just by their probability V {T ln{vo / v)) . 
For the exponential distribution we find that the "current- voltage" characteristics change from the Ohm 
low / ~ w at low V to the nonlinear regime / ~ y^/^o with a diverging differential resistance at lowest u, 
as shown on the Figure 123 

For the scaling distribution of the collective pinning we obtain a very slow decrease xZ^ ^ | lnti;|~" 
and / ^ |lnu|~". These results are very encouraging since they show the same functional forms as 
the formula H44(l for the collective pinning creep with a = 1/ ji = (2 — C)/x ^-^d the formula (|59|) with 
a = Ijx- Actually there is no discrepancy, even in powers, since Ea. H44|l was derived already for the 
periodic media, where C = 0. 

Certainly, the straightforward merging of results from collective and local pinning theories is spec- 
ulative and should be used only as a suggestion for more rigorous studies. Nevertheless the observed 
unification of different and and differently derived results of the collective pinning theory at 

one side, and essentially different view of the local pinning theory at another side, looks quite optimistic. 
We shall discuss some other aspects of this correspondence in the next Section. 

12 Interference of local and collective pinnings. 

We have already seen that the simplified but explicit approach of the local pinning provides clear inter- 
pretations for hypotheses of the collective pinning, particularly on the origin of metastable states. The 
same time, it raises challenges which have not been met yet. E.g. the following items of the above analysis 
are important for the theory of the w — / dependencies: 

1. A fraction of metastable branches terminate at the end points or relax fast at minimal points. 
There are those points which determine /(u) at high enough v but they are not accounted in the collec- 
tive pinning theory. 

2. A fraction of stronger metastable branches do not posses this instability which at first sight 
allows for the large v perturbative approach of the collective pinning theories. But it results in the even 
more obscure effect of generating sequences of dislocations or solitons. Now the v — f dependencies are 
determined by competing processes: the annihilation contrary to the aggregation. These processes are 
not accessible yet to existing theories except for the simple treatment of the local pinning which also 
needs to be further elaborated. Particularly demanding are studies of aggregation and annihilation of 
dislocations, their own pinning, etc. 

Consider more closely some other aspects of interference of collective and local pinning centers. Prag- 
matically, we shall concentrate on those which will be important in studies of the response functions 
in applications to susceptibility of density waves discussed in the next Section. The problem is more 
fundamental, being related to interference of different scales within the collective pinning picture. Here, 
a great simplification comes from the clear separation of both the length and the time scales between 
the local pinning and the collective pinning regimes. The slow evolution of the collective pinning allows 
us to consider it within a given distribution of local metastable states. Then the collective pinning is 
described by the same Hamiltonian H93fl where the generic pinning potential V{(p{'x.)) is substituted by 
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the local pinning energy V{ip{'K)) =^ Ea{i^{'X-) + Qx,). The collective pinning evolution will be described 
by the same equation (|39|l where the generic random force is substituted by Fa — E'^. within the phase 
language it reads 

i 

Without local effects, the sum over impurities would correspond to the random force density g(x, (^(x)) 
of the Eq. (|39|l . Here the major difference from the pure collective pinning is that now = i^^ is a 
two- valued function where the degree of freedom a = ±ora=: 1,2 allows for the thermodynamic or 
stochastic treatments, contrary to the frozen disorder x;. 

For the linear response, or for the motion with a small velocity, most of impurities are locally relaxed, 
and we shall neglect others for a moment. Then the random potential becomes single valued, a = 1 only, 
but its correlator i?, defined by the equation lO, changes in comparison with the generic form R ^ cos(ip). 
To make things simple, consider the limit of a very strong impurity potential V. Then the lowest energy 
state is V = V'i(^) ~ d — 2ttQ{6 — tt) and the pinning energy becomes W{ipi) which originates the 
discontinuous force. E.g. for the short range model of the Sec. I7.2.2l we have V{9) = — cos6') => V*{9) 
with 

V*{e) = W{Me))^EAl-siu^-^—^\ , F ^ ^ ^\nLgn{e - n) 



Remind that these energies and forces already take into account the elastic adaptation, so their 
correlator should be compared to the renormalized correlator of pinning forces rather than with the bare 
one . We see that the force becomes discontinuous at = tt which property is more general than our 
particular choice of the model. The bare smooth correlator (|10|l becomes singular, cusp-like, at (/s = 0: 

R= / —V*i9)V*{e + ^))^cos{^)^{\^\-n)cos^-2sm^-f 



Apparently it contains the nonanalytic therm ~ — which originates the kink ~ in the correlator 
of forces, 

- I (tt — \lp\) cos — — 2 sin — 
4 '^'^ 2 2 



i?^^ A((p) ^ - ( (tt - cos ^ - 2 sin ^ ) (122) 



We can already observe the apparent link to the cusp anomaly in the force-force correlator discussed 
in the Sec. 14. II as a clue to the threshold pinning force. In this way, the local pinning picture suggests a 
transparent view and straightforward interpretation for one of the most important results of the collective 
pinning theory obtained with the help of the ERG. 

The bare kink in the random force correlator, originated by the local pinning at L <C ip, provides the 
necessary boundary condition |76| for the kink formation, and hence the development of the threshold 
field, within the collective pinning domain L ^ Lp. Also we can get an interpretation that the kink is 
formed by choosing the lowest state every time when the retarded and advanced terms cross each other 
changing their character from stable to metastable and vice versa. 

As well as in the case of the collective pinning, the cusp is rounded, if T 7^ 0, for small Lp; here it 
happens at AE = Ft^ip < T when both levels a = 1, 2 become comparably populated. 

Consider now those pinning centers which are not in equilibrium; each of them provides a point 
pinning force which is not random: it is directed against the applied force /. These states are close to 
the degeneracy, 9 = (^(x) -I- Qxi « ±7r -|- 69 and their fraction ly = S9/2tt is small (as T/Ft^ for the linear 
response problem and as vt^^ for the slow sliding). Their concentration rt„e = i^rii = determines the 
distance L„e which is large in compare to the mean distance between impurities. Still it can be either 
larger or smaller than the pinning length Lp. If v is not too small, such that Li <C Lne ^ Lp, there 
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A((p)i 




Figure 26: Comparison of the force correlators A{(p) = —Ripip{(p) following from local and collective 
pinning, respectively. The upper curve shows the force correlator of the local pinning approach as given 
by equation (|122(l . The lower curve shows the fixed point function A*(</?) ~ (<p — tt)^ — 7r^/3 found from 
collective pinning Note that A{(p) is periodic: A{(p) — A{(p + 27rn), n integer. The scales on the 
vertical axis has been chosen differently for the two curves to allow for a better comparison. 

are many nonequilibrium impurities within the pinning volume, and their point forces add to the total 
restoring force fioc of the local pinning which we have been studying before. The collective pinning will 
react to the difference /* = / — fioc developing its own reaction /co/(/*)- Then for the linear responses to 
both forces, /{, = Xb^^^^ '^'^ fi^'i ^^'^ additivity of the inverse susceptibilities which will be an important 
element of applications (e.g. Eq. (|123|l '). In case of very low when Lp <C L^e, we face the picture of 
very distant point sources of nonrandom forces. The reaction of the pinned elastic media to the isolated 
point force is not quite known and we can only guess that they will still contribute additively to the 
average pinning force. 

13 Some applications to Density Waves. 

Experimental observations on sliding charge and spin density waves, are very rich and clean; most general 
effects are very stable and observed similarly in different materials [51IBl[7i ril3| . At high enough T, the 
collective pinning picture is well confirmed in general. A typical observation is the inverse relation between 
the critical field for the onset of sliding and the real part 5Re — e' of the dielectric susceptibility e ~ x- 
fc oc 1/e' ^ P E]. The collective pinning is always affected by the T dependence of the elastic 
moduli (e.g. via the order parameter vanishing near the transition temperature Tc of the CDW/SDW 
formation, or via screening of the long-range Coulomb interactions at low T), and these are readily 
identified experimentally |114| . The critical dependencies of parameters C and V on Tc — T are known 
microscopically, interestingly different for CDWs and SDWs, and their combination confirms in all cases 
the T dependence of the critical field fc- The sliding also demonstrates the expected saturation of the 
V — f dependence to the linear law at high v. The local pinning does not show up at these high T as it 
should be: the barriers U cannot be higher than Tc, which is the scale of Eg, so that the relaxation is too 
fast for any observations. 

The picture changes drastically (see jll5| for the modern review) at T low in comparison to Tc and to 
the activation energy A of normal carriers, the last one is important because of the Coulomb hardening of 
the elasticity C,, ~ exp(A/T). In addition to the usual threshold fc = fti, the v — f curve shows a sharp 
upturn at the "second threshold field" ft2 (for reviews on earlier observations see pi I117L [TTHI ). Even 
beyond details, the overall v — f curvature becomes opposite to the high T one and hence to expectations 
of the collective pinning theory. The dielectric susceptibility e ~ x starts to show ui and T dependencies, 
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with 9e showing the maximum as a function of to and 5Re showing a surprising sharp peak as a function 
of T . These changes can be related to the opposite T dependencies for strengths of coUcctive and 
local pinnings, with the last one playing the major role at low T and not very low to. (Even at low T 
the collective pinning reemerges at ultra low oj which appear in measurements of the time delayed heat 
response |12()| .'l 

Below, we shall apply the picture of metastable plastic deformations to interpret these observations. 
The same model will allow to describe the two remarkable features which became commonly observed 
in Charge and Spin Density Waves. There are both the anomalous peak of £'{T,uj = const) and the 
nonlinear current- voltage / — £" (that is i; — /) curve with the second threshold field in the sliding regime. 
Namely, the features of £ result from a competition of the local relaxation with the collective pinning 
affected by the freezing of the Coulomb screening. The upper critical field in I — £ curves is reached 
when the shortest life time configurations are accessed by the fast moving density wave. 

13.1 Nonlinear f(v). 

Apparently, the first critical field £ti can be only the threshold due to the collective pinning £ti ^ fc- It 
seems to be almost time-independent which requires for high barriers available only within the collective 
pinning regime. The slow creep, observed as a "broad band noise" at finite T, corresponds to a distribution 
of high barriers in accordance with the collective pinning picture. 

Contrarily, £t2 appears to be the high velocity limit of the pinning force via the energy dissipated by 
the moving density wave which we identify with the maximal force derived above for the local pinning: 
£t2 ^ fmax- With increasing v, we shall inevitably reach the local pinning regime with its lower barrier 
heights necessary to provide the condition wt ~ 1. The qualitative curve of the Figure ESI shows the 
necessary positive curvature and the approach to the almost vertical / — as it is observed sometimes. 

In recent experiments [1151 two distinct regimes have been established: the linear one / ~ f 

at small £ followed by the exponential growth / ~ exp[cnst £]. Next, very recent experiments |121j 
have shown that £t2 is a steep crossover, rather than a kink as it was supposed for a long time, and 
it is closer to our picture. Now it is possible to fit quantitatively, within the same set of parameters, 
the f — V dependence over several orders of magnitude of v encompassing three different regimes of the 
theory described in the Sec. IHlas la, lb and 2, that is f ^ v, Inv and / — fmax ~ —1/v. The Figure 
1271 shows such a fit |121| done by the general formula H106|) with functions U{6) and F{9) specified for 
the short-range model (#2 of the Section l7.2|l . Even the slowing down at high currents corresponds 
qualitatively to high velocity effects related to generation of dislocations. Sec. 1^1 

13.2 Low T, low Lo susceptibility peak. 

Sliding density waves are principally characterized by their giant dielectric susceptibility, e' = 5fte ^ 
10^ — 10^, corresponding to the low threshold field £t ~ fc- Remarkably, a sharp maximum of e'{uj, T) as 
a function of T has been observed in a wide range of density waves materials at low T and for very low 
frequencies w |119l 1122"] . With lo decreasing from, typically, \Q^Hz to lO'^i/z, the maximum height is 
growing while its position Tmax{^) is shifting towards low T as shown on the Figure OSI Importantly, the 
uprising parts T > Tmax of all plots for e'(u; = const, T) follow the same master curve £'(T) and differ 
only by the cutoff Tmaxif^) below which e' drops sharply, see the Figure!^ 

While some features of the T dependence (the rising high T slope) are quite specific to CDWs and 
SDWs, we suggest this example as showing simultaneously a combination of several important ingredients: 
sensitivity of the collective pinning to elastic parameters; separation of time scales between the two types 
of pinning, as well as their interference in observations. 



57 




Voltage [V] 



Figure 27: The single model fits (after Miyano and Ogawa) for the experimental v{f) (current-voltage) 
curve through three different regimes, la, lb, 3a of the Sec. 18.21 The middle part of the semi-log plot 
clearly demonstrates the f ^ Inv regime. The left part, if plotted in the normal scale, shows the linear 
/ ~ u law. These two regimes indicate the presence of bistable pinning centers of either intermediate 
or high strengths. The sharp upturn at higher / discriminates in favor of the unrestricted bistability 
(/ ~ fmax — const /v). The curvature sign changing at low v is an artifact of the logarithmic rescaling 
of the current axis. The decreasing growth rate at highest v agrees with expected effects of dislocations, 
Sec. El 

The electric polarization, being proportional to the average phase displacement Sip, creates the restor- 
ing force fpin, which may decay in time. The external electric field £ opposes fpm which has two well 
separated contributions fpin = fcoi + hoc- The inverse susceptibility can be defined as: 

-'^--^^-c-i+eroi (123) 

Keeping track only of the major dependencies on T, Ui, we can write with the help of \112\ 

e-' ^ const n^c-^^l^ + , 'T^ (124) 

l + (iwro)-iexp(-C/^/T) 

where the first and the second terms come from the collective and the local pinnings correspondingly. The 
activation law e~^l^ in the first term comes from the effect of the Coulomb hardening of the longitudinal 
elastic modulus which is special to quasi one-dimensional density waves. 

For the model with the exponential distribution of barriers of the Sec. 1111 we obtain 

= Aexp(-^) + Buj^l^° , ^ - n,^ B ~ (125) 

where the second coefficient B is a complex weak function of T^uj. In both cases, the function e'{T,uj = 
const) is not monotonic, it demonstrates a peak at T^axi'^) defined by the equation 

1 , EaUo Vq 

uj exp(--2 — + Inn,) (126) 

'''0 ^ max ^ max 
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Figure 28: The dielectric susceptibility of density waves from interference of the local and the collective 
pinnings. The plots show the T dependencies of e'[(jj,T) at various uj. Calculations have been done for 
the formula (|125|l with A/B = 9 and Uo = 0.5 (both Uq and T are in units of A). While units of e' and 
UJ arbitrary, their changes correspond to the experimental plots below: at each step, ui was rescaled by 
one order of magnitude resulting in the overall change of e' by one order. 

Numerical calculations for (|125(l , shown on the Figure EHl are in agreement with the results of jll9j 
shown on the Figure 1^ Concerning the frequency dependence, the function 3e^|^.(r = const, u) in 
l)124|l is monotonous, but the measured function '^£ow{T = const, oj) shows a maximum which is also in 
agreement with experiments. 

The interpretation of these results |29j is that different types of pinning compete to contribute to 
, so that the lowest e dominates and near Tmax the pinning force is minimal. Namely, at the higher 
T slope the T dependence arises from the dispersionless (very high barriers for collective metastable 
states) collective pinning affected by statically screened Coulomb interactions. At the lower T slope, the 
dispersion of e comes from the relaxation of local metastable states. With increasing T at given oj, the 
local states approach the thermal equilibrium; then fioc — > 0, hence e{T, lo) grows with increasing T until 
the collective pinning force becomes dominating. Recall another, more phenomenological interpretation 
|119| which suggested approaching a kind of a dipole glass transition^. 

We conclude that the crossover picture of the local and the collective types of pinning describes 
altogether the susceptibility anomaly, long time relaxation and nonlinear v — f characteristics. The peak 
of £{T) results from the competition of the local metastable plasticity with the collective pinning affected 
by freezing out of the Coulomb screening. In / — w curves, the upper threshold field /2 is reached when 
the metastable plastic deformations with shortest life time are accessed by the fast moving density wave. 

^Indeed, there is an interesting, while indirect link to physics of two-level systems in conventional polar glasses 1511 . Thus, 
^Toc interpreted as e^"^ = a* , where cr* = —iuie* are the effective complex "conductivity" cr* and the "dielectric 

susceptibility" e* per one effective two-level "dipole" with "polarization" it/-,r under the "random field" 8 — n and the 
"external field" SO. 
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Figure 29: Experimental data by Nad et al (119) . Labels at the top of each curve show the frequency, in 
Hz. 



14 Conclusions 

In this article we have suggested an overview of the present status of our understanding of pinning 
phenomena in various systems, stressing the unifying aspects. Both languages of the collective and the 
local pinning exploit intensively the concept of metastable states. 

The picture of collective pinning discussed in the first part of this review is to some extend worked out. 
The theory has been successfully applied already for the explanation of equilibrium phases, depinning and 
creep phenomena in vortex lattices, charge density waves and magnetic domain walls. The decisive step 
here were the understanding of thermodynamic and dynamic scaling behavior, as well as development 
of renormalization group methods which are capable to consider the effects of metastability emerging on 
large length scales even from weak pinning centers. The domain of application of this theory are large 
length and time scales, the critical field vicinity and the creep below it. The theoretical description of 
plasticity in disordered systems is still incomplete. 

The picture of the local pinning, within its domain of low temperatures and not too low frequencies 
or velocities, can be effectively used to explain experimental data: qualitatively and even quantitatively. 
The advantages come from the explicit treatment of metastable states, their creation and relaxation, 
their relation to plasticity and topological defects. It provides also a clue to the quantum creep showing 
that the tunneling repulsion of crossing branches destroys the pinning. 
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Table 2: Frequently appearing quantities 



Li.ci.±±ij J. \j y 


y 1I.I.U\J1 




field coupling to Vu 


A 


(|H3 


cutoff 


ao 




Burgers vector 


b 


SccEl 


stiffness constant 


c = (q|Ci)i/3 




elastic moduli 


C||, c± 


® 


curvature 


c 


m 


dimension of the elastic object 


D 




space dimension 


d 


© 


energy barrier on scale L 


Eb{L) 


m 


energy of the dislocation loop of the radius R 


Edl{R) 




soliton energy 


Es 


Figurfd SecETa 


energy difference between low lying states 


AE 




frequency dependent barrier 






energy branches (1: stable, 2: metastable, +: accending, -: 


El, 2, E± 




descending) 






free energy 


T{L,u) 
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APPENDICES. 



A Free energy fluctuations in D = 1 dimensions 

An illustrative example is given by a linear D = 1-dimcnsional object such as a magnetic flux line with 
the boundary conditions u(L) = u and tt(0) = 0. Changing u enforces the object to see another disorder 
environment. Using the transfer matrix technique, it can be shown that F{L,u) obeys the equation |45| 



dL 2C du^ C \ du 



+ Vr{L,u). (127) 



As a side remark we mention that if we read L, u and Tr as time, space and height coordinates, 
respectively, eq. (|127f) becomes the Kardar-Parisi-Zhang-equation, which describes the height profile of 
a growing surface under the random influx Vr{L^u) of particles |123| . The correlations of the restricted 
free energy show the following scaling behavior: 

(:f,^l.u) - = (^)" * (^) . (m 

For small values of the argument y = (u — u')/wr{L) of $ the difference of the free energies should not 
depend on L since the configurations dominating the free energy will be the same for most parts of the 
elastic object (apart from x close to L). This gives $(2/) '-^ y^^^- For large arguments y the elastic 

,2\l/2 



stiffness dominates over the disorder and {(Tr{L, u) — J-r{L, u')) j ^ ^" ^" — ^ as in pure systems. A 

numerical solution of H127|l shows that for large L and intermediate values of y the free energy !Fr{L,u) 
forms a rugged landscape as a function of u with typical valleys of width w(L) separated by hills of 
height Tp{L/Lp)^ |124) . The general picture of a rugged energy landscape as concluded from eq. H128|l 
is believed to hold also for higher dimensional elastic objects. 



B Strong pinning in D = 1— dimensional CDWs 

To give a specific example for the case of strong pinning, we consider a lattice model for a charge density 
wave with the Hamiltonian given by 

H= J2 ^Qj(w. -5I^^^o^(27r(u,-K,))) (129) 

<i,j> i 

where 27rKi = — QXj and is a random impurity positions at which Vi ^ 0. If we assume for simplicity 
Vi = V for all i and consider the limit V ^ oo, then Ui = Hi + ni with ni integer and the Hamiltonian 
can be rewritten as ^ 

i7 = ^ -Cij{ni-nj + n-i- (130) 

<^,3> 

The minimization of this Hamiltonian leads to a set of integers {?T-i,o} from which a well defined result 
for ground state and hence the roughness wr follows. The ground state consists of regions of constant rii 
separated by oriented domain walls at which rii changes by ±1. A very simple situation exists in c? = 1 
dimensions, where < z, j >= z, i + 1 and the ground state follows trivially as rii+i = + [k^+i — k^] . Here 
[...] denotes the Gauss bracket which replaces its argument by the closest integer. Thus the Ui undergo 
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a random walk and hence C = 1/2. For a more detailed discussion of the one-dimensional case see e.g. 

urn- 

The specific Transfer Matrix technics (Appendix^ in the dimension D = 1 allows for a more detailed 
description of interference between the pinning and the thermal motion |53| . Thus for temperatures high 
in compare to the characteristic elastic energy T ^ T* ^ Cuimp, and arbitrary with respect to V, the 
heat capacitance c{T) is 

{V/Tf , T > F 

In the same regime, the correlations of the order parameter cos(27rui) decay exponentially with the 
correlation length ^(T) such that = + (^'^ . Here ~ C/T is the correlation length of a pure 
system while the randomness contribution to is 

^R-n~^^p {h{V/T)lh{VlT)f 

where Im are the modified Bessel functions. 

For low temperatures with a constraint to the strong pinning regime T <^ T* <^ V the correlation 
function of displacements behaves as a kind of the Mott law 

/jn^N 1/4 / fT*\^^^\ 
{{\ut-Uj\)j)j^^\i- j\ni„,pi — \ exp I -CTist f — j 1 

This nontrivial T dependence appears because thermal jumps take place primarily within segments of an 
optimal spacing \i — j]^^^ ~ n~^p{T* /Ty^^ which is much larger than the typical one 

C Details on metastable branches. 

Here we give details to the results of Sec. [7|on the energy branches and their special points in the language 
of the CDWs and its phase. 

Consider first the termination points defined in H99(l : all quantities at this point will be labelled by 
the index e. Expanding at the vicinity of the end point 

^^tP^ + S^P; 6* = 6*6 + 69, 69 <0 (131) 

we find from and the solutions 

/ V" \ 9 ('9T/'M3/2 

5^3,2 - ± [^jpjm) , u = v^m'/' , K = -^jj^^ . H'j = + < (132) 

As a function of V, the coefficient Ve is singular at F = Vi when the end point emerges and aX V = V2 
when it annihilates with the next one, 9 ^ 9^: 

1. Consider the emergence of metastable branches when the points 9^, —9^ + 2tt split from the point 
TT. We find 

V ^Vi+6V >Vi : e,~TT^ {6Vf'^ , U - {6V)-^/^{-59f'^ , - {6VY'^ (133) 

2. Consider the crossover to the unrestricted bistability: V ^ V2 when the two sets of end points 
9e, 9e join together and with the point 9m of the minimal barrier. At V = V2 + 6V < V2 the degeneracy 
is lifted and the branch crossing point {9m,ipm) splits into two end points 

9e, 9e = 9mT 69e 69^ ^ U ~ (6Vfl^{-69fl'^ (134) 
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At V > V2, U{9) passes through the minimum Umin ~ {5VY/'^ &i 9 — Om- For both signs of 5V we 
can write an interpolation 

V ^¥2: U ^ {BiSe^ + B2f^^ ; Bi - SV^'"^ , B2 5V 

3. Consider the limit of very strong impurities which allows for an explicit treatment. The equation 
W(V') + V s\-n{'ip — 0) = at F 3> maxW^' has the following solutions: 

E+ : ^^9- V-^W'{9) , E+ « W{9) - {2V)-^W'^{9) 

E_ : ^f=i9-2T:-V-^W'{9-2Ti), E_f=iW{9 -271)- {2V)~^W''^ {9- 2ti) 

E3 : ^Ri9-TT + V~^W'{9 - tt) , £^3 w + W{9 - n) - V'^W''^{9) 

U = E:i-E+p^2V + W{9-T:)-W{9),W'{9m.-T^)^W{9ra) 

4. Fortunately, for a point impurity we can order the branches and simplify the energy a priori even 
at arbitrary V . In what follows, tt < 9 < 2tt, while W, W are functions of ijj at —2it < tp < 2tt. For 
each term, we determine its own function ij; = i)a{9) with a — { — , +, 3} = {1, 2, 3} (within the selected 
semiperiod of 9). 





= 1!)+ : 9 = 


. W 

i- arcsm 




= -.9 = ^) 


. W 

-f arcsin 


V' 


= ip3 : 9 = n -\ 


- tp ~ arcsin 


Ui9) 


= E,{M0))- 





W-V\ll--^ +v 



V 



V 



These expressions were the bases for our plots on Figures [181221 

5. Consider in more details the overshooting branches which appear due to special long range effects of 
dislocations (recall the end of the Section [7!^ and the Section l^^T)) . This is the regime of small Sip ^ ip — 2TT 
and 59 = 9 — 2tt. Here, V is close to its minimum 5V ~ b/2{Stp — 59)"^, b = V"{0); W is close to its 
maximum W27r = 2i?s, but the expansion is not analytical. We shall write it, according to (|113fl taken for 
d = 3, as SW ~ —4:/3a{—Sip)^/^ where a ~ CR for the dislocation loop of the radius R. The minimization 
of H over -0 gives S9 = 6ip + 2a/b{-6ipy/'^ = 0. 

At 69 < 0, there is one solution 50 = - {a/b+{{a/bf - 39^^^- it gives the branch E^ approaching 
the end of the period, 9 ^ 2tt — with some deficiency E+{2tt) < 0: SE^ = — (8/3)a^/6'^ corresponding 
to the retardation Sip = — (2a/5)^. 

At S9 > 0, there are two solutions Sip = — (a/b ± {{a/by — S9y^^) . Here the sign — corresponds to 
the barrier branch E^ , the sign + corresponds to the overshooting part E^ of the branch E+ . 

Entering the next circle 69 > 0, the energy E+ keeps increasing, passing through the energy 2Es 
at Sip= - {Sa/2bf , S9 = 3/4(a/6)2. Further on it crosses with the branch E^ to become metastable. 
Since then, the difference E^ — E'^ = U gives the relaxation barrier. Finally, the two solutions E+ and 
E3 collapse at the termination point 

S9e = {a/bf , SiPe = -{a/bf, SH^ = 20a*/35^ 

The above results give rise to the picture of the Figure!^ and related conclusions. 
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D Details on the kinetic equation. 



The kinetic equation is derived from the balance law for occupation numbers of branches +,— , see Eq, 
(TO : 

dTi+ dn 



W^n+ + W±n^ , — = W±n^ + W±n+ (135) 
itween the branches and the 

d _d -d ^ ■ de{t) 



dt " "+ + ' ' dt 

Here W± , W± are transition rates between the branches and the full time derivative is 



dt dt de' dt 



V 



(136) 



W^/W± = exp {AE/T) , W±, W^^ exp {-Eb/T), AE = E+ - E^ (137) 
The relaxation rate is 

+ ^ exp [^^-^ ^ j (138) 

where Tq^ is an attempt rate. For AE >> T, r ~ exp{U/T) with the activation energy U = E^ — E2- 
Notice that at the end points the metastable branch disappears and there must he t{9 6^) ^ 0. Still 
the expression pU4|l leaves us with small but finite value of t ^ tq even at U{6f.) — * 0. Hence it should 
be corrected to provide r at f7 ^ T; it happens via a dependence to(0) which plausible form is 

ro~ (0-0e)^ fc> 1. 



D.l Stationary motion. 



Consider a stationary process when the density wave moves at a constant phase velocity v = —ip = 9 = 
const, then dn/dt = 0. Now the solution of the Eq. Ij 104(1 is trivial, but the boundary conditions must be 
properly specified. Suppose first that there are no end points which is the case of very strong impurities 
V > V2, see the Figure [T^ Then for 6 approaching 2n, both branches + and — contribute to initial 
conditions for the branch + at = adding the pair of solitons at infinity, see the Figure^] Contrarily, 
there is no source for the branch — a.t 6 = 0. This conditions read 



de 



d_ 

'de 



(ri+ + n_)27r = , hence n(0) = n, 



eq 



The solution of l(TM|) . l(T^ is 

n = ne,(so)e''''^''' + / neq[si)e''^~'' dsi 
where s{9) is an effective t/r over the branch: 

s = s{6) = 



de 



So 



.s(0) = -s(27r) 



(139) 



(140) 



(141) 



(We shall keep the same notations for functions of 9 and of s = s(0).) At presence of end points ( right 
6e and left 27r — 6*5) there is only one branch of lowest energy Ei which survives beyond {2tt — 9^, 9e), 
that is 

< 6* < 271-6*6 : n^l; e^<9 <2n : n = -1 (142) 
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The contributions of these monostable regions to the pinning force are exactly compensated as it should 
be. Within the bistability region {2tt — 9^, 0e) the solution is 



= / ne<j(s - S2) exp(-S2)ds2 (143) 

JQ 

The substitution to the general expression for the force l| 105(1 yields (at presence of end points) 

dsAE / (1326^"' [no{s - S2) - no{s)] (144) 
-00 Jo 

Mostly we shall consider the case of low T when AE S> T in essential regions, then 

neqis) « -sgns, n - Ueq = 9(s)e~'' (145) 

(8 is the unit step function). The point of symmetrical population, n = 0, is shifted to s = In 2. At 
< s < In 2 there is an inverted population, 712 > ni as shown on the Fig. 1171 Finally the expression 
H144|l is simphfied to the form H106() . 

D.2 Various regimes for f{v). 

Small velocities for all cases. 

At fr^ <^ 1, s(9) is large almost everywhere, except for a vicinity of tt where the barrier activation 
energy takes its largest value {/(tt) = Utt; then we should use (|144|l . Namely, if for largest T(7r) = we 
have vTtt 1, then s ~ 1 already at 6* — tt fr^ ^ 1 so that at (6* — tt) ~ 1 we have s ^ 1, then the 
series in v is well convergent. In lowest order of vt we find 

/ « nn, r dsAE^-^ = nn, f ( ^) ' dO « vr^F^ (146) 

Jo ds^ 'J cosh2(Ai;/2r) V de J " " ' 

At low T the dependence of the expression under the last integral in H146|l is governed by the factor, see 
P38|l . exp [— ((i?3 — E2) — (i?2 — £^1)) /T]. It has a maximum exp(— J/tt/T)) at 9 = n and we arrive at 
the result ((Tn5|l . 

At higher velocities, still only the vicinities of the crossing point « tt is important, but we must 
take into account the reduction of the barrier with increasing 9: U = C/tt — ^^(6' — 7r)/27r. We obtain 

f dsexp(-s) _ r n,T^F^v at vt^F^T < 1 , , 

' J s + VT^F^/{2nT)^\ n,T\n{vT^F^/T) at t-r^F^/T » 1 ^ 



High velocities: restricted metastability. Let wt^ >> 1, then s ^ 1 only at 6* « 9^. The form 
H107|l is more appropriate for calculations. Since e^'' — 1 is small at s ^ 1 i.e. at almost all 9, then only 
a vicinity of 9e contribute, hence we can take F = F^ a.t 9 ~ 9^- We obtain 

/ = fmax - 2^n,Fe J^^ d9 [l - e-^C')) ; = 2nn,AEe (148) 

Remind that at 9e the activation U = E3 — E2 ^ vanishes while the force is finite F^ = dAE/d9 7^ 0. 
We find 

fmax I Fe,TLi60y 
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where 69y is defined by the condition 

s{se,) ^ ^se-'^+^eM-^K) - i; ve= (f)'^'^ - (i49) 

and we arrive at the result H110() . vaUd in this simple form at Ueexp(— Vg/T) <^ v ^ Ve- 
High velocities: unrestricted metastability. 

The calculations are similar to the above case of the restricted metastability and we shall skip equiv- 
alent steps. The difference is that now there is a high velocity range v 3> Um = maxr~^ where the 1/v 
expansion is valid: 



/ = 2n, }^AE{2n) --J _(Ai?(2^) - AE{9))j 



(150) 



D.3 Linear response. 

Consider as a perturbation in the kinetic equation and expand as n = nc,q{9) + 5n{9,t): 



^ + i") 5„ + e^n,q = 0- f = n, I d9F{9)5n{9) (151) 



In the Fourier representation we have 



d_ 



U = - : -7Y / d9—j^—neg-- ■ (153) 



[-l+{^^r)-') A d9 d9 [^-l+i^u;r)-') 



Xc^ = T7- = -I . -I — , (154) 

which confirms (|112|1 . 

References 

[1] Monceau, P., 1985, in " Electronic properties of quasi one- dimensional compounds. '\ P. Monceau ed., 
Reidcl Publ. Co. 

[2] Griiner, G. 1994, Density Waves in Solids", Addison- Wesley Publ. Company, Reading. 

[3] Gor'kov, L., and Griiner, G., 1989, "Charge Density Waves in Solids'", Elsevier Sci. Publ., Amster- 
dam, 6, 365. 

[4] Dumas, J., and Schlenker, C, 1993, Int. Journ. Mod. Phys. B 7, 4045. 

[5] Proceedings of the International Workshop on Electronic Crystals, June 1993, ed. by Brazovskii, S., 
Monceau, P., 1993, J. Phys. France IV, CoUoquc C2 3. 

[6] Proceedings of the International Workshop on Electronic Crystals, 1999, ed. by Brazovskii, S., Mon- 
ceau, P., 1999, J. Phys. France IV, 9. 



69 



[7] Proceedings of the International Workshop on Eleetronic Crystals, September 2002, ed. by Bra- 
zovskii, S., Kirova, N., and Monceau, P., 2002, J. Phys. France IV, 12. 

X. Zhu, P. B. Littlewood, and A. J. Millis, 1994, Phys. Rev. B 50, 4600; and rfs. theirin. 

Giamarchi, T., and Orignac, E., in New Theoretical Approaches to Strongly Correlated Systems, ed. 
by A.M. Tsvelik, Kluwer 2001. 

Blatter, G., Feigel'man, M. V., Geshkenbein, V. B., Larkin, A. I. and Vinokur, V. M., 1994, Rev. 
Mod. Phys., 66, 1125. 

Nattermann, T., and Scheidl, S., 2000, Adv. in Physics, 49, 607-704. 

Lyuksyutov, I., Naumovets, V. and Pokrovsky, V., 1992, Two-dimensional Crystals, Academic Press, 
Boston. 

see, e.g.. Young, P.A., (ed.), 1999, Spin Glasses and Random Fields, World Scientific, Singapoore. 
Fisher, D.S., 1998, Phys. Rep. 301, 113. 
Kardar, M., 1998, Phys. Rep. 301, 85. 

Haasen, H.P., 1958, in Berichte der Cottinger Akademie der Wissenschaften. 
Larkin, A.I., 1970, Sov. Phys. JETP 31, 784. 

Larkin, A.I., and Ovchinnikov, Yu.N., 1979, J. Low Temp. Phys. 34, 409. 
Fukuyama, H., 1976, J. Phys. Soc. Jpn. 41, 513, 542. 
Fukuyama, H., and Lee, P.A., 1978, Phys. Rev. B 17, 535. 
Lee, P.A., and Rice, T.M., 1979, Phys. Rev. B 19, 3970. 
Nattermann, T., 1990, Phys. Rev. Lett. 64, 2454. 
Korshunov, S.E., 1993, Phys. Rev. B 48, 3969. 

Giamarchi, T., and Lc Doussal, P., 1994, Phys. Rev. Lett. 72, 1530; 1995, Phys. Rev. B 52, 1242. 
Fisher, D.S., 1985, Phys. Rev. B 31, 7233. 
Fisher, D.S., 1986, Phys. Rev. Lett. 56, 1964. 
Caroli, C., in |7j, p. 269, and rfs. therein. 
Larkin, A.I., 1994, Sov. Phys. JETP 105, 1793. 
Larkin, A.I., and Brazovskii, S., 1995, Sol. State Commun. 93, 275. 
Brazovskii, S., in |113j . p. 467. 

Brazovskii, S., and Larkin, A.I., in [Sj and cond-mat / 99 11100 



Nattermann, T., 2000, in Spin Glasses and Random Fields, ed. by Young. A. P.. World Scientific, 
Singapoore. 



70 



[33] Nattcrmann, T., 1985, Z. Physik, B54,247 (1984), Phys. stat. sol. (B), 132, 125. 
[34] Emig, T., and Nattcrmann, T., 1998, Phys. Rev. Lett. 81, 1469. 
[35] Grinstcin, G., and Ma, S.-K., 1982, Phys. Rev. Lett. 49, 685. 

[36] Landau, L.D., and Lifshitz, L.F., 1989, Theoretical Physics, Vol. 7, Akademie-Vcrlag Berhn. 
[37] Efetov, K.B., and Larkhi, A.I., 1977, Sov. Phys. JETP. 45, 1236. 
[38] Andrei, E.Y., et aL, 1988, Phys. Rev. Letters 60, 2765. 
[39] WiUet, R.L., et al., 1989, Phys. Rev. B 38, R 7881. 

[40] Giamarchi, T., 2002, in Proceedings of the NATO ASI "Field Theory of Strongly Correlated Fermions 
and Bosons in Low- Dimensional Disordered Systems", Kluwer. 

[41] Imry, Y. and Ma, S.K., 1975, Phys. Rev. Lett. 35, 1399. 

[42] Alava, M.J., et al., 2001, in "Phase Transitions and Critical Phenomena", ed. by C. Domb and J.L. 
Lebowitz, Vol. 18, p. 145, Academic Press. 

[43] Villain, J., and Semeria, B., 1983, J. Phys. (Paris), 44, L 889. 

[44] Villain, J., and Fernandez, J.F., 1984, Z. Physik B 54, 139. 

[45] Huse, D.A., Henley, C.L., and Fisher, D.S., 1985, Phys. Rev. Lett. 55, 2924. 

[46] Nattcrmann, T., 1987, Europhysics Lett. 4, 1241; Nattcrmann, T., 1988, Phys. Rev. Lett 60, 2701. 

[47] Mezard, M., 1990, J. Phys. (Paris) 51, 1831. 

[48] Fisher, D.S., and Huse, D.A., 1991, Phys. Rev. B43, 10728. 

[49] Schulz, U., Villain, J., Brczin, E., and Orland, H. 1988, J. Stat. Phys., 51, 1. 

[50] Korshunov, S.E., 2001, Phys. Rev. B 63, 174514. 

[51] Anderson, P.W., Halperin, B.I., and Varma, C.M., 1972, Phil. Mag. 8 (25), 1. 
[52] Schehr, G., Giamarchi, T., and Le Doussal, P., 2003, Phys. Rev. Lett. 91, 117002. 
[53] Brazovskii, S., 1979, Sov. Phys.: JETP 49, 504. 

[54] Patashinsky, A.Z. and Pokrovsky, V., 1979, "Fluctuation Theory of Phase Transitions", Pergamon 
Press Oxford. 

[55] Amit, D.J., 1984, "Field Theory, the Renormalization Group, and Critical Phenomena", World 
Scientific Singapore. 

[56] Balents, L. and Fisher, D.S., 1993, Phys. Rev. B 48, 5949. 

[57] Balents, L., Bouchaud, J. P., and Mezard, M., 1996, J. Phys. I (Paris) 6 1007. 

[58] Chauve, P., Le Doussal, P., and Wiese, J.-K., 2001, Phys. Rev. Lett. 86, 1785. 

[59] Lassig, M., 1998, Phys. Rev. Lett. 80, 889. 



71 



[60] For a review see also Halpin-Healy, T.. and Zhang. Y.C., 1995, Physics Reports 254, 215. 

[61] Emig, T., Bogner, S., and Nattermann, T., 1999, Phys. Rev. Lett. 83, 400; S. Bogner, T. Emig and 
T. Nattermann, 2001, Phys. Rev. B 63 , 174501. 

[62] Kicrfcld, J., Hwa, T., and Nattermann, T., 1997, Phys. Rev. B. 55, 626. 

[63] Carpentier, D., Le Doussal, P., and Giamarchi, T., 1996, Europhys. Lett. 35, 15193. 

[64] Fisher, D.S., 1997, Phys. Rev. Lett. 78, 1964. 

[65] Kierfeld, J., 1998, Physica C 300, 171. 

[66] Kierfeld, J., and Vinokur, V., 2000, Phys. Rev. B 61, R14928. 

[67] Feigel'man, M., 1980, Sov. Phys. 52, 555. 

[68] Glatz, A., and Nattermann, T., 2002, Rev. Lett. 88, 256401. 

[69] Littlewood, P.B., and Rammal, R., 1988, Phys. Rev. B 38, 2675; Myers, C.R., and Sethna, J.P., 
1993, Phys. Rev. B 47, 11194. 

[70] Drosscl, B., and Kardar, M., 1995, Phys. Rev. E, 52, 4841. 

[71] Mikheev, L.V., Drossel, B., and Kardar, M., 1995, Phys. Rev. Lett., 75, 1170. 

[72] loffe, L. B., and Vinokur, V. M.,1987, J. Phys. C 20, 6149. 

[73] Nattermann, T., 1987, Europhys. Lett, 4, 1241. 

[74] Radzihovsky, L., 1998, Talk E37 8 given at the March Meeting of the American Physical Society. 

[75] Chauve, P., Giamarchi, T., and Le Doussal, P., 1998, Europhys. Lett. 44, 110. 

[76] Narayan, O., and Fisher, D.S., 1992, Phys. Rev. Lett., 68, 3615; Phys. Rev. B, 46, 11520. 

[77] Nattermann, T., Stepanow, S., Tang, L.-H., and Leschhorn, H., 1992, J. Phys. II France 2, 565 
(1483); 1997, Ann. Phys. (Leipzig) 6, 1. 

[78] Narayan, O., and Fisher, D.S., 1993, Phys. Rev. B 48, 7030. 

[79] Ertas, D., and Kardar, M., 1994, Phys. Rev. E 49, R2532. 

[80] Vinokur, V.M., and Nattermann, T., 1997, Phys. Rev. Lett. 79, 3471. 

[81] Glatz, A., and Nattermann, T., 2003, |cond-mat/030405'0 



[82] Fisher, D.S., 1983, Phys. Rev. Lett. 50, 1486. 
[83] Middleton, A.A., 1991, Phys. Rev. B 45, 9465. 



[84] Roters, E.L., Usabel, K.D., and Liibeck, S. cond-mat/0207137 and references therein. 
[85] Glatz, A., Nattermann, T., and Pokrovsky, V., 2003, Phys. Rev. Lett. 90, 047201. 
Kohandel, M. and Kardar, M. 1999, Phys. Rev. B 59, 9637. 



72 



[87] KohandcL M. and Kardar, M., 2000, Phys. Rev. B 61, 11729. 

[88] Zeng, C, Leath, P. L., and Fisher, D. S., 1999, Phys. Rev. Lett., 82, 1935. 

[89] Fricdel, J., 1964, "Dislocations", Addison- Wesley, Reading, Mass.. 

[90] Nabarro, F., 1967, "Theory of Crystal Dislocations", Oxford University, New York. 

[91] Dumas, J., and Feinberg, D., 1986, Europhys. Lett. 2, 555. 

[92] Feinberg, D., and Friedel, J., 1988, J. Phys. France 49, 485. 

[93] Brazovskii S., in O , p. 425. 

[94] Brazovskii, S., and Matveenko, S., 1991, Sov. Phys.: JETP 72, 492; ibid, 72, 860. 
[95] Nad', F.Ya., in p.l91. 

[96] Gor'kov, L.P., 1983, JETP Letters 38, 87 and in gj. 

[97] Ong, N.P. and Maki, K., 1985, Phys. Rev. B 32, 6582. 

[98] GiU, J.C., 1986, J. Phys C 19, 6589; in Ref. 0, p. 165; in Ref. 0, p. 165. 

[99] Langer, J.S., and Ambegoakar, V., 1967, Phys.Rev. 164, 498. 

[100] Ivlev, B.I., Kopnin, N.B., 1984, Advances in Physics 33, 47. 

[101] Varoquaux, E., and Avenel, O., 1994, Physica B 197, 306. 

[102] Itkis, M.E., Nad', F.Ya., Monceau, P., and Renard, M., in [S], p.l75. 

[103] Lemay, S.G., et al, 1998, Phys. Rev. B 57, 12781 and in 0, 1993, p. 149; and rfs. therein 

[104] Requardt, H. et al, 1998, Phys. Rev. Lett., 80, 5631. 

[105] Brazovskii, S., and Kirova, N., in p. 139. 

[106] Coppersmith S.N., Phys. Rev. Lett., 1990, 65, 1044. 

[107] Rideau, D., 2001, et al, Europhysics Letters 56, 289. 

[108] Kobelkov A., Gleisberg F. and Wonneberger W., 1999, Eur. Phys. J. B, 9, 21. 

[109] Brazovskii S. and Kirova N. in [7], p. 173. 

[110] Zaitsev-Zotov, S.V., 1993, Phys. Rev. Lett., 71, 605. 

[Ill] Ji-Min Duan, Phys. Rev. B 48, 4860 ; Maki, K., 1995, Phys. Lett. A, 202, 313; Matsukawa, H., 
Miyake, H., Yumoto, M., Fukuyama H. and Kawaguchi T., in |^, p. 161. 

[112] Caldeira, A.O. and Leggett, A.J., 1981, Phys. Rev. Lett. 46, 211. 

[113] Physics and Chemistry of Low Dimensional Conductors, 1996, Schlenker C. et al eds.. Proceedings 
of the NATO- as: series B: Physics, vol. 354, Plenum Press. 



73 



[114] Maki K. and Virosztck A., Phys. Rev. B, 1989, 39, 9640; ibid, 1990, 42, 655; Tomic, S., Cooper, 
J.R., Kang, W., Jrome, D. and Maki, K., J.Physique I, 1991, 1, 1603. 

[115] Thome, R.E., et al in 0, p. 291. 

[116] Lcmay, S.G., Thornc, R.E., Li Y., and Brock J.D., Phys. Rev. Lett. 83, 2793 (1999). 
[117] Maeda, A., Notomi, M., Kim, T.W., Uchinokura, K., 1990, Phys. Rev. B 42, 3290. 
[118] Itkis, M., Nad, F.Ya., Monccau, P., 1990, J. Phys.: Condensed matter 2, 8327. 
[119] Nad, F.Ya., and Monceau, P., 1993, Sohd State Commun. 87, 13; and in 0, p. 343. 
[120] Biljakovic, K., 1998, in and in Physica B, 1993, 244, 167. 

[121] Miyano, K. and Ogawa, N., 2001, in [TJ; Ogawa, N., and Miyano, K., 2001, Phys. Rev. Lett., 2001, 
87, 256401; and Phys. Rev. B, submitted. 

[122] Lasjaunias, J.C., Biljakovic, K., Nad, F.Ya., Monceau, P., Bechgaard, K., 1994, Phys. Rev. Lett. 72, 
1283. 

[123] Kardar, M., Parisi, C, and Zhang, Y-C, 1986, Phys. Rev. Lett. 56, 889. 

[124] Kardar, M., 1987, J. AppL Phys. 61, 3601; Kardar, M., and Zhang, Y.C., 1987, Phys. Rev. Lett. 
58, 2087. 

[125] Glatz, A., Kumar, S., and Li, M.S., 2001, Phys. Rev. B 64, 184301. 



74 



